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Abstract 

In this article we review the present status of the spin foam formulation of non- 
perturbative (background independent) quantum gravity. The article is divided in 
two parts. In the first part we present a general introduction to the main ideas 
emphasizing their motivations from various perspectives. Riemannian 3-dimensional 
gravity is used as a simple example to illustrate conceptual issues and the main goals 
of the approach. The main features of the various existing models for 4-dimensional 
gravity are also presented here. We conclude with a discussion of important questions 
to be addressed in four dimensions (gauge invariance, discretization independence, 
etc.). 

In the second part we concentrate on the definition of the Barrett-Crane model. 
We present the main results obtained in this framework from a critical perspective. 
Finally we review the combinatorial formulation of spin foam models based on the 
dual group field theory technology. We present the Barrett-Crane model in this 
framework and review the finiteness results obtained for both its Riemannian as well 
as its Lorentzian variants. 
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1 Introduction 



Quantum gravity, the theory expected to reconcile the principles of quantum mechanics 
and general relativity, remains a major challenge in theoretical physics (for a review of the 
history of quantum gravity see pQ). The main lesson of general relativity is that, unlike 
in any other interaction, space-time geometry is fully dynamical. This special feature of 
gravity precludes the possibility of representing fields on a fixed background geometry 
and severely constrains the applicability of standard techniques that are successful in the 
description of other interactions. Although the necessity of a background independent 
formulation of quantum gravity is widely recognized, there is a current debate about the 
means by which it should be implemented. In particular, it is not clear whether the non 
renormalizability of perturbative quantum gravity should be interpreted as an indication of 
the inconsistency of general relativity at high energies, the inconsistency of the background- 
dependent framework applied to gravity, or a combination of both. 

According to the (background dependent) perspective of standard QFT 2j, non renor- 
malizability signals the inconsistency of the theory at high energies to be corrected by a 
more fundamental theory in the UV regime. A classical example of this is Fermi's four- 
fermion theory as an effective description of the weak interaction. According to this view 
different approaches to quantum gravity have been defined in terms of modifications of 
general relativity based on supersymmetry, higher dimensions, strings, etc. The finiteness 
properties of the perturbative expansions (which are background dependent from the on- 
set) are improved in these theories; however, the definition of a background independent 
quantization of such modifications remains open. 

The approach of non perturbative quantum gravity is based on a different interpre- 
tation of the infinities in perturbative quantum gravity: it is precisely the perturbative 
(background dependent) techniques which are inconsistent with the fundamental nature of 
gravity. This view is strongly suggested by the predictions of the background independent 
canonical quantization of general relativity defined by loop quantum gravity. Loop quan- 
tum gravity (LQG) is a non perturbative formulation of quantum gravity based on the 
connection formulation of general relativity (for an reviews on the subject see jSl IU El E] ) • 
A great deal of progress has been made within the theory. At the mathematical level, the 
main achievement is the rigorous definition of the Hilbert space of quantum geometry, the 
regularization of geometric operators and the rigorous definition of the quantum Hamilto- 
nian constraint (defining the quantum dynamics). States of quantum geometry are given 
by polymer-like excitations supported on graphs (spin network states). From the physical 
viewpoint its main prediction is the discreteness of geometry at the Planck scale. This 
provides a clear-cut understanding of the problem of UV divergences in perturbative gen- 
eral relativity: at the Planck scale the classical notion of space and time simply ceases to 
exist; therefore, it is the assumption of a fixed smooth background geometry (typically flat 
space-time) in perturbation theory that becomes inconsistent at high energies. The theory 
successfully incorporates interactions between quantum geometry and quantum matter in 
a way that is completely free of divergences [Zj. The quantum nature of space appears as 
a physical regulator for the other interactions. 
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Dynamics is governed by the quantum Hamiltonian constraint. Even when this op- 
erator is rigorously defined jH] it is technically difficult to characterize its solution space. 
This is partly because the 3 + 1-decomposition of space-time (necessary in the canonical 
formulation) breaks the manifest 4-diffeomorphism invariance of the theory making awk- 
ward the analysis of dynamics. The situation is somewhat analogous to that in standard 
quantum field theory. In the Hamiltonian formulation of standard quantum field theory 
manifest Lorentz invariance is lost due to a particular choice of time slicing of Minkowski 
space-time. The formalism is certainly Lorentz invariant, but one has to work harder to 
show it explicitly. Manifest Lorentz invariance can be kept only in the Lagrangian (path- 
integral) quantization making the (formal) path integral a powerful device for analyzing 
relativistic dynamics. 

Consequently, there has been growing interest in trying to define dynamics in loop 
quantum gravity from a 4-dimensional covariant perspective. This has given rise to the 
so-called spin foam approach to quantum gravity. Its main idea is the construction of a 
rigorous definition of the path integral for gravity based on the deep insights obtained in 
the canonical framework of loop quantum gravity. In turn, the path integral provides a 
device to explicitly solve the dynamics: path-integral transition amplitudes can be shown 
to correspond to solutions of the quantum Hamiltonian constraint. 

The underlying discreteness discovered in loop quantum gravity is crucial: in spin foam 
models the formal Misner-Hawking functional integral for gravity is replaced by a sum 
over combinatorial objects given by foam-like configurations (spin foams). A spin foam 
represents a possible history of the gravitational field and can be interpreted as a set of 
transitions through different quantum states of space. Boundary data in the path integral 
are given by the polymer-like excitations (spin network states) representing 3-geometry 
states in loop quantum gravity. General covariance implies the absence of a meaningful 
notion of time and transition amplitudes are to be interpreted as defining the physical 
scalar product. 

While the construction can be explicitly carried out in three dimensions there are 
additional technical difficulties in four dimensions. Various models have been proposed. A 
natural question is whether the infinite sums over geometries defining transition amplitudes 
would converge. In fact, there is no UV problem due to the fundamental discreteness and 
potential divergences are associated to the IR regime. There are recent results in the 
context of the Barrett-Crane model showing that amplitudes are well defined when the 
topology of the histories is restricted in a certain way. 

The aim of this article is to provide a comprehensive review of the progress that has 
been achieved in the spin foam approach over the last few years and provide as well a 
self contained introduction for the interested reader that is not familiar with the subject. 
The article is divided into two fundamental parts. In the first part we present a general 
introduction to the subject including a brief summary of LQG in Section |21 We introduce 
the spin foam formulation from different perspectives in Section El In Section 0] we present 
a simple example of spin foam model: Riemannian 3-dimensional gravity. We use this 
example as the basic tool to introduce the main ideas and to illustrate various conceptual 
issues. We review the different proposed models for 4-dimensional quantum gravity in Sec- 
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tion Finally, in Section |H] we conclude the first part by analyzing the various conceptual 
issues that arise in the approach. The first part is a general introduction to the formalism; 
it is self contained and could be read independently. 

One of the simplest and most studied spin foam model for 4-dimensional gravity is 
the Barrett-Crane model [HJ US]. The main purpose of the second part is to present a 
critical survey of the different results that have been obtained in this framework and its 
combinatorial generalizations fTJ EI] based on the dual group field theory (GFT) 
formulation. In Section [7| we present a systematic derivation of the Barrett-Crane model 
from the £pin(4) Plebanski's formulation. This derivation follows an alternative path from 
that of here we emphasize the connection to a simplicial action. 

Spin foams can be thought of as Feynman diagrams. In fact a wide class of spin foam 
models can be derived from the perturbative (Feynman) expansion of certain dual group 
field theories (GFT) ^SJ IT7j . A brief review of the main ideas involved is presented in 
Section |S1 We conclude the second part by studying the GFT formulation of the Barrett- 
Crane model for both Riemannian and Lorentzian geometry. The definition of the actual 
models and the sketch of the corresponding finiteness proofs [T%irH?ll2U] are given in Section 

El 

Part I 

Main ideas 

2 Loop Quantum Gravity and Quantum Geometry 

Loop quantum gravity is a rigorous realization of the quantization program established in 
the 60's by Dirac, Wheeler, De-Witt, among others (for recent reviews see jHllUE]). The 
technical difficulties of Wheeler's 'geometrodynamics' are circumvent by the use of con- 
nection variables instead of metrics fH\ 1221 12H] • At the kinematical level, the formulation 
is similar to that of standard gauge theories. The fundamental difference is however the 
absence of any non-dynamical background field in the theory. 

The configuration variable is an £[/(2)-connection A l a on a 3-manifold £ representing 
space. The canonical momenta are given by the densitized triad Ef. The latter encode 
the (fully dynamical) Riemannian geometry of E and are the analog of the 'electric fields' 
of Yang-Mills theory. 

In addition to diffeomorphisms there is the local £77(2) gauge freedom that rotates the 
triad and transforms the connection in the usual way. According to Dirac, gauge freedoms 
result in constraints among the phase space variables which conversely are the generating 
functionals of infinitesimal gauge transformations. In terms of connection variables the 
constraints are 

Q l = V a E? = 0, C a = E b k E b k a = 0, S = (? jh E?E b j F abk + --- = 0, (1) 
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where T> a is the co variant derivative and Fb a is the curvature of A l a . Qi is the famil- 
iar Gauss constraint — analogous to the Gauss law of electromagnetism — generating in- 
finitesimal SU(2) gauge transformations, C a is the vector constraint generating space- 
diffeomorphism, and S is the scalar constraint generating 'time' reparameterization (there 
is an additional term that we have omitted for simplicity). 

Loop quantum gravity is defined using Dirac quantization. One first represents (0) as 
operators in an auxiliary Hilbert space 7i and then solves the constraint equations 

g^ = o, 4^ = o, sv = o. (2) 

The Hilbert space of solutions is the so-called physical Hilbert space H p h ys - In a generally 
covariant system quantum dynamics is fully governed by constraint equations. In the case 
of loop quantum gravity they represent quantum Einstein's equations. 

States in the auxiliary Hilbert space are represented by wave functionals of the connec- 
tion ^(A) which are square integrable with respect to a natural diffeomorphism invariant 
measure, the Ashtekar-Lewandowski measure (we denote it C?[A] where A is the space 
of (generalized) connections). This space can be decomposed into a direct sum of orthog- 
onal subspaces Ti = @ 7 7Y 7 labeled by a graph 7 in E. The fundamental excitations are 
given by the holonomy he(A) G SU{2) along a path £ in E: 

hi{A) = V exp / A. (3) 

Elements of 7i 7 are given by functions 

¥ /i7 (A) = f{heM), M^)> • • • . M^)). ( 4 ) 

where hi is the holonomy along the links I G 7 and / : SU(2) n — > C is (Haar measure) 
square integrable. They are called cylindrical functions and represent a dense set in TC 
denoted Cyl. 

Gauge transformations generated by the Gauss constraint act non-trivially at the end- 
points of the holonomy, i.e., at nodes of graphs. The Gauss constraint (in (JIJ) is solved by 
looking at SU(2) gauge invariant functionals of the connection (C 2 [A\/Q). The fundamen- 
tal gauge invariant quantity is given by the holonomy around closed loops. An orthonormal 
basis of the kernel of the Gauss constraint is defined by the so called spin network states 
^7,Of},W}(^) E3 IlHl HZI • Spin-networks 1 are defined by a graph 7 in E, a collection of 
spins {je} — unitary irreducible representations of SU(2) — associated with links I G 7 and 
a collection of SU{2) intertwiners {i n } associated to nodes n G 7 (see Figure [T]). The spin- 
network gauge invariant wave functional \l/ 7 i{j } M (A) is constructed by first associating 
an SU(2) matrix in the j^-representation to the holonomies hg(A) corresponding to the 

1 Spin- networks were introduced by Penrose [28) in a attempt to define 3-dimensional Euclidean quantum 
geometry from the combinatorics of angular momentum in QM. Independently they have been used in 
lattice gauge theory 00] as a natural basis for gauge invariant functions on the lattice. For an account 
of their applications in various contexts see |31) . 
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link £, and then contracting the representation matrices at nodes with the corresponding 
intertwiners t n , namely 

%.m,M^) = ll n hMA)}, (5) 

nG7 €£7 

where ji[h e (A)} denotes the corresponding ^-representation matrix evaluated at corre- 
sponding link holonomy and the matrix index contraction is left implicit. 



5_ 
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Figure 1: Spin-network state: At 3-valent nodes the intertwiner is uniquely specified by the 
corresponding spins. At 4 or higher valent nodes an intertwiner has to be specified. Choos- 
ing an intertwiner corresponds to decompose the n-valent node in terms of 3-valent ones 
adding new virtual links (dashed lines) and their corresponding spins. This is illustrated 
explicitly in the figure for the two 4-valent nodes. 

The solution of the vector constraint is more subtle [21]. One uses group averaging 
techniques and the existence of the diffeomorphism invariant measure. The fact that 
zero lies in the continuous spectrum of the diffeomorphism constraint implies solutions to 
correspond to generalized states. These are not in 7i but are elements of the topological 
dual Cyl* 2 . However, the intuitive idea is quite simple: solutions to the vector constraint 
are given by equivalence classes of spin-network states up to diffeomorphism. Two spin- 
network states are considered equivalent if their underlying graphs can be deformed into 
each other by the action of a diffeomorphism. 

This can be regarded as an indication that the smooth spin-network category could 
be replaced by something which is more combinatorial in nature so that diffeomorphism 
invariance becomes a derived property of the classical limit. LQG has been modified 
along these lines by replacing the smooth manifold structure of the standard theory by 
the weaker concept of piecewise linear manifold [32] ■ In this context, graphs defining 
spin-network states can be completely characterized using the combinatorics of cellular 
decompositions of space. Only a discrete analog of the diffeomorphism symmetry survives 
which can be dealt with in a fully combinatorial manner. We will take this point of view 
when we introduce the notion of spin foam in the following section. 

2 According to the triple Cyl dHc Cyl*. 



8 



Quantum geometry 

The generalized states described above solve all of the constraints $1} but the scalar con- 
straint. They are regarded as quantum states of the Riemannian geometry on S. They 
define the kinematical sector of the theory known as quantum geometry. 

Geometric operators acting on spin network states can be defined in terms of the fun- 
damental triad operators E®. The simplest of such operators is the area of a surface S 
classically given by 

A S {E) = [ dx 2 ^Ti[n a n h E^E b ] (6) 
Js 

where n is a co normal. The geometric operator A${E) can be rigorously defined by its 
action on spin network states [33| l3~H 185] . The area operator gives a clear geometrical 
interpretation to spin-network states: the fundamental 1-dimensional excitations defining 
a spin-network state can be thought of as quantized 'flux lines' of area. More precisely, 
if the surface 5* C £ is punctured by a spin-network link carrying a spin j, this state 
is an eigenstate of As(E) with eigenvalue proportional to j(j + 1). In the generic 
sector — where no node lies on the surface — the spectrum takes the simple form 

a s ({j}) = 8tu4 VWi + 1)> (7) 

i 

where i labels punctures and i is the Imirzi parameter 3 - a s{{j}) is the sum of single 
puncture contributions. The general form of the spectrum including the cases where nodes 
lie on S has been computed in closed form 

The spectrum of the volume operator is also discrete [33J El EH EH] • If we define the 
volume operator V a (E) of a 3-dimensional region a C S then non vanishing eigenstates are 
given by spin- networks containing n-valent nodes in a for n > 3. Volume is concentrated 
in nodes. 

Quantum dynamics 

In contrast to the Gauss and vector constraints, the scalar constraint does not have a simple 
geometrical meaning. This makes its quantization more involved. Regularization choices 
have to be made and the result is not unique. After Thiemann's first rigorous quantization 
[iU] other well defined possibilities have been found [JTJ H21 EH] • This ambiguity affects 
dynamics governed by 

= 0. (8) 

The difficulty in dealing with the scalar constraint is not surprising. The vector 
constraint — generating space diffeomorphisms — and the scalar constraint — generating time 
reparameterizations — arise from the underlying 4-diffeomorphism invariance of gravity. In 

3 The Imirzi parameter i is a free parameter in the theory. This ambiguity is purely quantum mechanical 
(it disappears in the classical limit). It has to be fixed in terms of physical predictions. The computation 
of black hole entropy in LQG fixes the value of i (see |37|'l. 



9 



the canonical formulation the 3 + 1 splitting breaks the manifest 4-dimensional symmetry. 
The price paid is the complexity of the time re-parameterization constraint S. The sit- 
uation is somewhat reminiscent of that in standard quantum field theory where manifest 
Lorentz invariance is lost in the Hamiltonian formulation . 

From this perspective, there has been growing interest in approaching the problem 
of dynamics by defining a covariant formulation of quantum gravity. The idea is that 
(as in the QFT case) one can keep manifest 4-dimensional covariance in the path integral 
formulation. The spin foam approach is an attempt to define the path integral quantization 
of gravity using what we have learn from LQG. 

In standard quantum mechanics path integrals provide the solution of dynamics as 
a device to compute the time evolution operator. Similarly, in the generally covariant 
context it provides a tool to find solutions to the constraint equations (this has been 
emphasized formally in various places: in the case of gravity see for example jH], for a 
detailed discussion of this in the context of quantum mechanics see (35] )• We will come 
back to this issue later. 

Let us finish by stating some properties of S that do not depend on the ambiguities 
mentioned above. One is the discovery that smooth loop states naturally solve the scalar 
constraint operator [IdI I4T] . This set of states is clearly to small to represent the physical 
Hilbert space (e.g., they span a zero volume sector). However, this implies that S acts only 
on spin network nodes. Its action modifies spin networks at nodes by creating new links 
according to Figure El 5 - This is crucial in the construction of the spin foam approach of 
the next section. 



k k 




i i 



Figure 2: A typical transition generated by the action of the scalar constraint 



3 Spin Foams and the path integral for gravity 

The possibility of defining quantum gravity using Feynman's path integral approach has 
been considered long ago by Misner and later extensively studied by Hawking, Hartle and 
others [1EJI3H]- Given a 4-manifold M. with boundaries Si and E 2 , and denoting by G the 

4 There is however an additional complication here: the canonical constraint algebra does not reproduce 
the 4-diffeomorphism Lie algebra. This complicates the geometrical meaning of S. 

5 This is not the case in all the available definitions of the scalar constraints as for example the one 
defined in g2J|13]. 
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space of metrics on M., the transition amplitude between \[q a b}) on Si and \[q' ab }) on S 2 is 
formally 



where the integration on the right is performed over all space-time metrics up to 4- 
diffeomorphisms [g] G G/ixff(M) with fixed boundary values up to 3-diffeomorphisms [q a b], 
Kb], respectively 

There are various difficulties associated with (JHJ). Technically there is the problem of 
defining the functional integration over [g] on the RHS. This is partially because of the 
difficulties in defining infinite dimensional functional integration beyond the perturbative 
framework. In addition, there is the issue of having to deal with the space G / Diff{M), i.e., 
how to characterize the diffeomorphism invariant information in the metric. This gauge 
problem (3-diffeomorphisms) is also present in the definition of the boundary data. There 
is no well defined notion of kinematical state \[q a b\) as the notion of kinematical Hilbert 
space in standard metric variables has never been defined. 

We can be more optimistic in the framework of loop quantum gravity. The notion of 
quantum state of 3-geometry is rigorously defined in terms of spin-network states. They 
carry the diff-invariant information of the Riemannian structure of S. In addition, and 
very importantly, these states are intrinsically discrete (colored graphs on E) suggesting 
a possible solution to the functional measure problem, i.e., the possibility of constructing 
a notion of Feynman 'path integral' in a combinatorial manner involving sums over spin 
network world sheets amplitudes. Heuristically, '4-geometries' are to be represented by 
'histories' of quantum states of 3-geometries or spin network states. These 'histories' 
involve a series of transitions between spin network states (Figure EJ), and define a foam- 
like structure (a '2-graph' or 2-complex) whose components inherit the spin representations 
from the underlying spin networks. These spin network world sheets are the so-called spin 
foams. 

The precise definition of spin foams was introduced by Baez in [T3] emphasizing their 
role as morphisms in the category defined by spin networks 6 . A spin foam T : s — > s', 
representing a transition from the spin-network s = (7, {je}, {t n }) into s' = (7', {je}, {v}), 
is defined by a 2-complex 7 J bordered by the graphs of 7 and 7' respectively, a collection 
of spins {jf} associated with faces / G J and a collection of intertwiners {i e } associated to 
edges e G J . Both spins and intertwiners of exterior faces and edges match the boundary 
values defined by the spin networks s and s' respectively. Spin foams T : s — > s' and 
T' : s' — > s" can be composed into TT' : s — > s" by gluing together the two corresponding 
2-complexes at s'. A spin foam model is an assignment of amplitudes A[F] which is 
consistent with this composition rule in the sense that 



6 The role of category theory for quantum gravity had been emphasized by Crane in |5UI 1511 152] . 
7 In most of the paper we use the concept of piecewise linear 2-complexes as in Jl] ; in Section [S] we 
shall study a formulation of spin foam in terms of certain combinatorial 2-complexes. 




(9) 



A\TT'\ = A[F]A[F']. 



(10) 
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Transition amplitudes between spin network states are defined by 

(s,s\ hys = Yl A m (ii) 

where the notation anticipates the interpretation of such amplitudes as defining the physical 
scalar product. The domain of the previous sum is left unspecified at this stage. We shall 
discuss this question further in Section El This last equation is the spin foam counterpart 
of equation 0. This definition remains formal until we specify what the set of allowed 
spin foams in the sum are and define the corresponding amplitudes. 




Figure 3: A typical path in a path integral version of loop quantum gravity is given 
by a series of transitions through different spin-network states representing a state of 3- 
geometries. Nodes and links in the spin network evolve into 1-dimensional edges and faces. 
New links are created and spins are reassigned at vertexes (emphasized on the right). The 
'topological' structure is provided by the underlying 2-complex while the geometric degrees 
of freedom are encoded in the labeling of its elements with irreducible representations and 
intertwiners. 

In standard quantum mechanics the path integral is used to compute the matrix ele- 
ments of the evolution operator U(t). It provides in this way the solution for dynamics 
since for any kinematical state \1/ the state U(t)^ is a solution to Schrodinger's equation. 
Analogously, in a generally covariant theory the path integral provides a device for con- 
structing solutions to the quantum constraints. Transition amplitudes represent the matrix 
elements of the so-called generalized 'projection' operator P (Sections 13.11 and 16. 3|) such 
that P\l/ is a physical state for any kinematical state As in the case of the vector 
constraint the solutions of the scalar constraint correspond to distributional states (zero 
is in the continuum part of its spectrum). Therefore, TC p hy S is not a proper subspace of H 
and the operator P is not a projector (P 2 is ill defined) 8 . In Section 0] we give an explicit 
example of this construction. 

The background-independent character of spin foams is manifest. The 2-complex can be 
thought of as representing 'space-time' while the boundary graphs as representing 'space'. 

In the notation of the previous section states in 7i p hy S are elements of Cyl* . 
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They do not carry any geometrical information in contrast with the standard concept of a 
lattice. Geometry is encoded in the spin labelings which represent the degrees of freedom 
of the gravitational field. 



3.1 Spin foams and the projection operator into Tt phys 

Spin foams naturally arise in the formal definition of the exponentiation of the scalar 
constraint as studied by Reisenberger and Rovelli in [23] and Rovelli The basic idea 
consists of constructing the 'projection' operator P providing a definition of the formal 
expression 

P=H S(S(x)) = j V[N\ e iS[N \ (12) 

where S[N] = j dx 3 N(x)S(x), with N(x) the lapse function. P defines the physical scalar 
product according to 

(s,s') phya = (sP,s'), (13) 

where the RHS is defined using the kinematical scalar product. Reisenberger and Rovelli 
make progress toward a definition of ()12j) by constructing a truncated version Pa (where A 
can be regarded as an infrared cutoff). One of the main ingredients is Rovelli's definition 
of a diffeomorphism invariant measure T>[N] generalizing the techniques of [55J. 
The starting point is the expansion of the exponential in ()12j) in powers 

\N(x)\<A \ n=0 ■ / 

The construction works for a generic form of quantum scalar constraint as long as it acts 
locally on spin network nodes both creating and destroying links (this local action generates 
a vertex of the type emphasized in Figure |2]). The action of S[N] depends on the value 
of the lapse at nodes. Integration over the lapse can be performed and the final result is 
given by a power series in the cutoff A, namely 



DC 



n=0 jF„:s— >s' 

= E l — E UMpv,^), (is) 

n=0 ' Fn.s^s' v 

where T n : s — > s' are spin foams generated by n actions of the scalar constraint, i.e., spin 
foams with n vertices. The spin foam amplitude A[jF n ] factorizes in a product of vertex 
contributions A v (p v , i v ) depending of the spins j v and i v neighboring faces and edges. The 
spin foam shown in Figure El corresponds in this context to two actions of S and would 
contribute to the amplitude in the order A 2 . 
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Physical observables can be constructed out of kinematical operators using P. If O kin 
represents an operator commuting with all but the scalar constraint then O phys = PO kin P 
defines a physical observable. Its expectation value is 

I ^ \ \ S PO kin P,s) \SP\O kin P\, S) /-,r>\ 

(sP, S) A^oo {sP A , S) 

where the limit of the ratio of truncated quantities is expected to converge for suitable 
operators O kin . Issues of convergence have not been studied and they would clearly be 
regularization dependent. 



3.2 Spin foams from lattice gravity 

Spin foam models naturally arise in lattice-discretizations of the path integral of grav- 
ity and generally covariant gauge theories. This was originally studied by Reisenberger 
[25] . The space-time manifold is replaced by a lattice given by a cellular complex. The 
discretization allows for the definition of the functional measure reducing the number of 
degrees of freedom to finitely many. The formulation is similar to that of standard lattice 
gauge theory. However, the nature of this truncation is fundamentally different: back- 
ground independence implies that it cannot be simply interpreted as a UV regulator (we 
will be more explicit in the sequel). 

We present a brief outline of the formulation for details see [221 I2E1 121] • Start from the 
action of gravity in some first order formulation (S(e, A)). The formal path integral takes 
the form 

Z = J V[e] V[A] e iS{e ' A) = J V[A] e iS ^ A \ (17) 

where in the second line we have formally integrated over the tetrad e obtaining an effective 
action S cB 9 . From this point on, the derivation is analogous to that of generally covariant 
gauge theory. The next step is to define the previous equation on a 'lattice'. 

As for Wilson's action for standard lattice gauge theory the relevant structure for the 
discretization is a 2-complex J . We assume the 2-complex to be defined in terms of the dual 
2-skeleton J7~a of a simplicial complex A. Denoting the edges e G J7a and the plaquettes 
or faces / G J7a one discretizes the connection by assigning a group element g e to edges 

A - {g e }. 

The Haar measure on the group is used to represent the connection integration: 

V[A] - H dg e . 

eeJ A 

9 This is a simplifying assumption in the derivation. One could put the full action in the lattice [SEj 
and then integrate over the discrete e to obtain the discretized version of the quantity on the right of (fTTjl . 
This is what we will do in Section QJ 
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The action of gravity depends on the connection A only through the curvature F(A) so 
that upon discretization the action is expressed as a function of the holonomy around faces 
gf corresponding to the product of the g e 's which we denote gf — gi 1} ■ ■ ■ g ( e n) - 

HA) -> {g f }. 

In this way, S eS (A) — > S eS ({gf}). Thus the lattice path integral becomes: 

Z= I II dg e exp[tS cS ({g f })}. (18) 

Reisenberger assumes that S cB ({gf}) is local in the following sense: the amplitude of 
any piece of the 2-complex obtained as its intersection with a ball depends only on the 
value of the connection on the corresponding boundary. Degrees of freedom communicate 
through the lattice connection on the boundary One can compute amplitudes of pieces of 
J7~a (at fixed boundary data) and then obtain the full J7a amplitude by gluing the pieces 
together and integrating out the mutual boundary connections along common boundaries. 
The boundary of a portion of J7a is a graph. The boundary value is an assignment of 
group elements to its links. The amplitude is a function of the boundary connection, 
i.e., an element of Cyl. In the case of a cellular 2-complex there is a maximal splitting 
corresponding to cutting out a neighborhood around each vertex. If the discretization is 
based on the dual of a triangulation these elementary building blocks are all alike and 
denoted atoms. Such an atom in four dimensions is represented in Figure EJ 

2 



Figure 4: A fundamental atom is defined by the intersection of a dual vertex in J7a (corre- 
sponding to a 4-simplex in A) with a 3-sphere. The thick lines represent the internal edges 
while the thin lines the intersections of the internal faces with the boundary. They define 
the boundary graph denoted 75 below. One of the faces has been emphasized. 

The atom amplitude depends on the boundary data given by the value of the holonomies 
on the ten links of the pentagonal boundary graph 75 shown in the figure. This amplitude 
can be represented by a function 

V(aij) for ctij G G and i ^ j — 1, . . . , 5 (19) 
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where a^- represents the boundary lattice connection along the link ij in Figure 0] Gauge 
invariance (V(ay) = V (diOiijgJ 1 )) implies that the function can be spanned in terms of 
spin networks functions ^-y B ,{pij}Mi}( a v) based on the pentagonal graph 75, namely 

V K) = EE^^^^) *7 6 ,{p^,mK-) (2°) 

Pij 

where V({pij}, {ti}) is the atom amplitude in 'momentum' space depending on ten spins 
Pij labeling the faces and five intertwiners tj labeling the edges. Gluing the atoms together 
the integral over common boundaries is replaced by the sum over common values of spin 
labels and intertwiners 10 . The total amplitude becomes 

z im= e e n Mpf,^) n v({p„},m), (21) 

C/l/H{p/}c e :{eH{ te }/ejA vej A 

where C e : {e} — > {i e } denotes the assignment of intertwiners to edges, Cf : {p/} — > {/} 
the assignment of spins pf to faces, and Af(pf,if) is the face amplitude arising in the 
integration over the lattice connection. The lattice definition of the path integral for 
gravity and covariant gauge theories becomes a discrete sum of spin foam amplitudes! 

3.3 Spin foams for gravity from BF theory 

The integration over the tetrad we formally performed in (|17j) is not always possible in 
practice. There is however a type of generally covariant theory for which the analog 
integration is trivial. This is the case of a class of theories called BF theory. General 
relativity in three dimensions is of this type. Consequently, BF theory can be quantized 
along the lines of the previous section. BF spin foam amplitudes are simply given by 
certain invariants in the representation theory of the gauge group. We shall study in some 
detail the case of 3-dimensional Riemannian quantum gravity in the next section. 

The relevance of BF theory is its close relation to general relativity in four dimensions. 
In fact, general relativity can be described by certain BF theory action plus Lagrange 
multiplier terms imposing certain algebraic constraints on the fields This is the 

starting point for the definition of several of the models we will present in this article: 
a spin foam model for gravity can be defined by imposing restrictions on the spin foams 
that enter in the partition function of the BF theory. These restrictions are essentially the 
translation into representation theory of the constraints that reduce BF theory to general 
relativity. 

10 This is a consequence of the basic orthogonality of unitary representations, namely 

/S 
j[9] a /3 k [9hndg = ^j^S^S^. 
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3.4 Spin foams as Feynman diagrams 



As already pointed out in ^3] spin foams can be interpreted in close analogy to Feynman 
diagrams. Standard Feynman graphs are generalized to 2-complexes and the labeling of 
propagators by momenta to the assignment of spins to faces. Finally, momentum conser- 
vation at vertices in standard Feynmanology is now represented by spin-conservation at 
edges, ensured by the assignment of the corresponding intertwiners. In spin foam models 
the non-trivial content of amplitudes is contained in the vertex amplitude as pointed out 
in Sections 13 . II and 13.21 which in the language of Feynman diagrams can be interpreted as 
an interaction. We shall see that this analogy is indeed realized in the formulation of spin 
foam models in terms of a group field theory (GFT) [TBI I17j. 



4 Spin foams for 3- dimensional gravity 

Three dimensional gravity is an example of BF theory for which the spin foam approach 
can be implemented in a rather simple way. Despite its simplicity the theory allows for the 
study of many of the conceptual issues to be addressed in four dimensions. In addition, 
as we mentioned in Section E31 spin foams for BF theory are the basic building block of 
4-dimensional gravity models. For a beautiful presentation of BF theory and its relation 
to spin foams see [SH]- For simplicity we study the Riemannian theory; the Lorentzian 
generalization of the results of this section have been studied in [60 . 

4.1 The classical theory 

Riemannian gravity in 3 dimensions is a theory with no local degrees of freedom, i.e., no 
gravitons. Its action is given by 

S(e,A) = J Tr(e A F(A)), (22) 

M 

where the field A is an S"?7(2)-connection and the triad e is a Lie algebra valued 1-form. 
The local symmetries of the action are SU(2) gauge transformations 

Se = [e, to] , 5 A = cIauj, (23) 

where u is a su(2)-valued 0-form, and 'topological' gauge transformation 

Se = d A r], 5 A = 0, (24) 

where d A denotes the covariant exterior derivative and rj is a su(2)-valued 0-form. The 
first invariance is manifest from the form of the action, while the second is a consequence 
of the Bianchi identity, d A F = 0. The gauge freedom is so big that the theory has only 
global degrees of freedom. This can be checked directly by writing the equations of motion 

F(A) = d A A = 0, d A e = 0. (25) 
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The first implies that the connection is flat which in turn means that it is locally gauge 
(A = d A u). The solutions of the second equation are also locally gauge as any closed form 
is locally exact (e = d A rj) 11 . This very simple theory can be quantized in a direct manner. 



4.2 Canonical quantization 

Assuming M = Exl where £ is a Riemann surface representing space, the phase space 
of the theory is parameterized by the spatial connection A (for simplicity we use the same 
notation as for the space-time connection) and its conjugate momentum E. The constraints 
that result from the gauge freedoms an d are 



The first is the familiar Gauss constraint (in the notation of equation ([!])), and we call the 
second the curvature constraint. There are 6 independent constraints for the 6 components 
of the connection, i.e., no local degrees of freedom. The kinematical setting is analogous to 
that of 4-dimensional gravity and the quantum theory is defined along similar lines. The 
kinematical Hilbert space (quantum geometry) is spanned by SU(2) spin-network states 
on E which automatically solve the Gauss constraint. "Dynamics" is governed by the 
curvature constraint. The physical Hilbert space is obtained by restricting the connection 
to be flat and the physical scalar product is defined by a natural measure in the space of flat 
connections |59j . The distributional character of the solutions of the curvature constraint 
is manifest here. Different spin network states are physically equivalent when they differ 
by a null state (states with vanishing physical scalar product with all spin network states). 
This happens when the spin networks are related by certain skein relations. One can 
reconstruct H p h ys directly from the skein relations which in turn can be found by studying 
the covariant spin foam formulation of the theory. 

4.3 Spin foam quantization of 3d gravity 

Here we apply the general framework of Section 13.21 This has been studied by Iwasaki in 
|611 162j . The partition function, Z, is formally given by 12 



where for the moment we assume Ai to be a compact and orientable. Integrating over the 
e field in (J2*T|) we obtain 



11 One can easily check that the infinitesimal diffeomorphism gauge action Se = C v e, and 5A = £ V A, 
where C v is the Lie derivative in the v direction, is a combination of ill' Ml and 124(1 for ui — v a A a and 
-q = v a e a , respectively, acting on the space of solutions, i.e. when (|2"Kll holds. 

12 We are dealing with Riemannian 3-dimensional gravity. This should not be confused with the approach 
of Euclidean quantum gravity formally obtained by a Wick rotation of Lorentzian gravity. Notice the 
imaginary unit in front of the action. The theory of Riemannian quantum gravity should be regarded as 
a toy model with no obvious connection to the Lorentzian sector. 



V a E1 = 0, F(A) = 0. 



(26) 




(27) 
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The partition function Z corresponds to the 'volume' of the space of flat connections on 
M. 

In order to give a meaning to the formal expressions above, we replace the 3-dimensional 
manifold M. with an arbitrary cellular decomposition A. We also need the notion of the 
associated dual 2-complex of A denoted by J7a- The dual 2-complex J A is a combinatorial 
object defined by a set of vertices v G J A (dual to 3-cells in A) edges e G J A (dual to 
2-cells in A) and faces / G J A (dual to 1-cells in A). 

The fields e and A have support on these discrete structures. The sw(2)-valued 1-form 
field e is represented by the assignment of an e G su(2) to each 1-cell in A. The connection 
field A is represented by the assignment of group elements g e G SU(2) to each edge in J A . 

The partition function is defined by 



Z(A) = 




where dej is the regular Lebesgue measure on R 3 , dg e is the Haar measure on SU(2), and 
U f denotes the holonomy around faces, i.e., Uf = g\ . . . g* for iV being the number of edges 
bounding the corresponding face. Since Uf G SU (2) we can write it as Uf = u°f 1 + Ff 
where «' 6 C and Ff G su(2). Ff is interpreted as the discrete curvature around the 
face /. Clearly Tr[e fUf] = Tr[e/F/]. An arbitrary orientation is assigned to faces when 
computing Uf. We use the fact that faces in J A are in one-to-one correspondence with 
1-cells in A and label e/ with a face subindex. 
Integrating over e/, we obtain 

Z(A)= j J] dg e J] <7f), (30) 

eeJ A feJ A 

where 5 corresponds to the delta distribution defined on C 2 (SU{2)). Notice that the 
previous equation corresponds to the discrete version of equation (J2EJ)- 

The integration over the discrete connection (fj dg e ) can be performed expanding first 
the delta function in the previous equation using the Peter- Weyl decomposition [63J 

5(g) = A ^[j(9)}, (31) 

iGirrep(SC/(2)) 

where Aj = 2j + 1 denotes the dimension of the unitary representation j, and j(g) is 
the corresponding representation matrix. Using equation ()31|). the partition function (J3TIJ) 
becomes 

^( A ) = E / II d 9e II Tr IM ■ ■ ■ 9e)} , (32) 

c-.{j}^{f} J eeJ A feJ A 

where the sum is over coloring of faces in the notation of (|2ip. 

Going from equation (j2U|) to we have replaced the continuous integration over the 
e's by the sum over representations of SU(2). Roughly speaking, the degrees of freedom 
of e are now encoded in the representation being summed over in (J32J) . 



19 



Now it remains to integrate over the lattice connection {g e }- If an edge e e J~ A bounds 
n faces there are n traces of the form Tr[jy(- • • g e ■ ■ ■ )] in (}3*2*)) containing g e in the argument. 
The relevant formula is 

PL ■= [ dg ./.(//) his) ® • • ■ o^fo) = E ^-*. ^-.v ( 33 ) 

where is the projector onto Invfjx g) j 2 £5 • • • <8> j n ]. On the RHS we have chosen 
an orthonormal basis of invariant vectors (intertwiners) to express the projector. Notice 
that the assignment of intertwiners to edges is a consequence of the integration over the 
connection. This is not a particularity of this example but rather a general property of 
local spin foams as pointed out in Section l3~2l Finally fpTTTf) can be written as a sum over 
spin foam amplitudes 

z(A) = £ En n Mivjv), (34) 

where is given by the appropriate trace of the intertwiners l v corresponding 

to the edges bounded by the vertex and j v are the corresponding representations. This 
amplitude is given in general by an SU (2) 3iVj-symbol corresponding to the flat evaluation 
of the spin network defined by the intersection of the corresponding vertex with a 2-sphere. 
When A is a simplicial complex all the edges in J7a are 3-valent and vertices are 4-valent 
(one such vertex is emphasized in Figure E3 the intersection with the surrounding S* 2 is 
shown in dotted lines). Consequently, the vertex amplitude is given by the contraction of 
the corresponding four 3-valent intertwiners, i.e., a 6j'-symbol. In that case the partition 
function takes the familiar Ponzano-Regge [64] form 




were the sum over intertwiners disappears since dim(Inv[7i ® j% <E> J3]) = 1 for 577(2) and 
there is only one term in (J33|) . Ponzano and Regge originally defined the amplitude (|35|) 
from the study of the asymptotic properties of the 6j-symbol. 

4.3.1 Discretization independence 

A crucial property of the partition function (and transition amplitudes in general) is that 
it does not depend on the discretization A. Given two different cellular decompositions A 
and A' (not necessarily simplicial) 

r- no Z(A) = r^Z(A'), (36) 
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where no is the number of O-simplexes in A (hence the number of bubbles in J7a), and 
r = + 1) is clearly divergent which makes discretization independence a formal 

statement without a suitable regularization. 

The sum over spins in (J35|) is typically divergent, as indicated by the previous equation. 
Divergences occur due to infinite volume factors corresponding to the topological gauge 
freedom (j2U)(see [65J) 13 . The factor r in represents such volume factor. It can also be 
interpreted as a 5(0) coming from the existence of a redundant delta function in (J3Uj) . One 
can partially gauge fix this freedom at the level of the discretization. This has the effect 
of eliminating bubbles from the 2-complex. 

In the case of simply connected S the gauge fixing is complete. One can eliminate 
bubbles and compute finite transition amplitudes. The result is equivalent to the physical 
scalar product defined in the canonical picture in terms of the delta measure 14 . 

In the case of gravity with cosmological constant the state-sum generalizes to the 
Turaev-Viro model [66j defined in terms of SU q (2) with q n = 1 where the representa- 
tions are finitely many. Heuristically, the presence of the cosmological constant introduces 
a physical infrared cutoff. Equation (JHUj) has been proved in this case for the case of sim- 
plicial decompositions in [SH], see also [HUE!]- The generalization for arbitrary cellular 
decomposition was obtained in |69j . 

13 For simplicity we concentrate on the Abelian case G — U(l). The analysis can be extended to the 
non-Abelian case. Writing g € U(l) as g — e l9 the analog of the gravity simplicial action is 

S(A,{e f },{0e})= e /^/({^})' ( 37 ) 

where Ff({9 e }) = J2 e ef® e - Gauge transformations corresponding to l|23|) act at the end points of edges 
e G c7a by the action of group elements {/?} in the following way 

B f B h 

9 e ^ d e + (3 S - (3 t , (38) 

where the sub-index s (respectively t) labels the source vertex (respectively target vertex) according to the 
orientation of the edge. The gauge invariance of the simplicial action is manifest. The gauge transformation 
corresponding to l|24|) acts on vertices of the triangulation A and is given by 

B f Bf + rjs-rjt, 

Oe^de- (39) 

According to the discrete analog of Stokes theorem 

f&Bubble 

which implies the invariance of the action under the transformation above. The divergence of the corre- 
sponding spin foam amplitudes is due to this last freedom. Alternatively one can understand it from the 
fact that Stokes theorem implies a redundant delta function in H30J) per bubble in J7a ■ 

14 If M. = S 2 x [0, 1] one can construct a cellular decomposition interpolating any two graphs on the 
boundaries without having internal bubbles and hence no divergences. 
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4.3.2 Transition amplitudes 



Transition amplitudes can be denned along similar lines using a manifold with boundaries. 
Given A, J7a then defines graphs on the boundaries. Consequently, spin foams induce 
spin networks on the boundaries. The amplitudes have to be modified concerning the 
boundaries to have the correct composition property (fTU|) . This is achieved by changing 
the face amplitude from (Aj ) to (A^) 1 ' 2 on external faces. 

The crucial property of this spin foam model is that the amplitudes are independent of 
the chosen cellular decomposition |?)7tl69j. This allows for computing transition amplitudes 
between any spin network states s = (7, {j}, {1}) and s' = (7, {j'}, {V}) according to the 
following rules 15 : 

• Given M. = E x [0, 1] (piecewise linear) and spin network states s = (7, {j}, {t}) and 
s ' = (7) {i'}> i 1 '}) 011 the boundaries — for 7 and 7' piecewise linear graphs in E— 
choose any cellular decomposition A such that the dual 2-complex J7a is bordered 
by the corresponding graphs 7 and 7' respectively (existence can be shown easily). 

• Compute the transition amplitude between s and s f by summing over all spin foam 
amplitudes (rescaled as in (jSHJ)) for the spin foams F : s — ► s' defined on the 2- 
complex J7a- 

4.3.3 The generalized projector 

We can compute the transition amplitudes between any element of the kinematical Hilbert 
space H 16 . Transition amplitudes define the physical scalar product by reproducing the 
skein relations of the canonical analysis. We can construct the physical Hilbert space by 
considering equivalence classes under states with zero transition amplitude with all the 
elements of H, i.e., null states. 



15 Here we are ignoring various technical issues in order to emphasize the relevant ideas. The most 
delicate is that of the divergences due to gauge factors mentioned above. For a more rigorous treatment 
see [TH] , 

16 The sense in which this is achieved should be apparent from our previous definition of transition 
amplitudes. For a rigorous statement see |70|. 




Figure 5: Elementary spin foams used to prove skein relations. 
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Here we explicitly construct a few examples of null states. For any contractible Wilson 
loop in the j representation the state 




V = (2j +l)s- {J ) ®s = a 0, (40) 



for any spin network state s, has vanishing transition amplitude with any element of 7i. 
This can be easily checked by using the rules stated above and the portion of spin foam 
illustrated in Figure El to show that the two terms in the previous equation have the same 
transition amplitude (with opposite sign) for any spin-network state in 7i. Using the second 
elementary spin foam in Figure one can similarly show that 



phys 



(41) 



or the re-coupling identity using the elementary spin foam on the right of Figure 03 

1 




^V2s + lV2e+ 1 



i j s 
k I e 



e =0, (42) 



where the quantity in brackets represents an SU(2) 6j-symbol. All skein relations can 
be found in this way. The transition amplitudes imply the skein relations that define the 
physical Hilbert space! The spin foam quantization is equivalent to the canonical one. 

4.3.4 The continuum limit 

Recently Zapata [7U] formalized the idea of a continuum spin foam description of 3- 
dimensional gravity using projective techniques inspired by those utilized in the canonical 
picture [21]. The heuristic idea is that due to the discretization invariance one can define 
the model in an 'infinitely' refined cellular decomposition that contains any possible spin 
network state on the boundary (this intuition is implicit in our rules for computing transi- 
tion amplitudes above). Zapata concentrates on the case with non- vanishing cosmological 
constant and constructs the continuum extension of the Turaev-Viro model. 



4.4 Conclusion 

We have illustrated the general notion of the spin foam quantization in the simple case of 
3 dimensional Riemannian gravity (for the generalization to the Lorentzian case see [oTIj). 
The main goal of the approach is to provide a definition of the physical Hilbert space. The 
example of this section sets the guiding principles of what one would like to realize in four 
dimensions. However, as should be expected, there are various new issues that make the 
task by far more involved. 
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5 Spin foams for 4- dimensional quantum gravity 



In this section we briefly describe the various spin foam models for quantum gravity in the 
literature. 



5.1 The Reisenberger model 

According to Plebanski |HH] the action of self dual Riemannian gravity can be written as 
a constrained SU(2) BF theory 



S(B, A) = J Tr [B A F(A)} — ipi 



1'/ A IP - -5 ij B k A B k 
3 



(43) 



where variations with respect to the symmetric (Lagrange multiplier) tensor ip^ imposes 
the constraints 

Q ij = B i A B j - -5 ij B k AB k = 0. (44) 
3 

When B is non degenerate the constraints are satisfied if and only if B % = ±(e°Ae l +|e J jfc e J ' A 
e h ) which reduces the previous action to that of self-dual general relativity. Reisenberger 
studied the simplicial discretization of this action in jSHj as a preliminary step toward the 
definition of the corresponding spin foam model. The consistency of the simplicial action 
is argued by showing that the simplicial theory converges to the continuum formulation 
when the triangulation is refined: both the action and its variations (equations of motion) 
converge to those of the continuum theory. 

In reference jSZ] Reisenberger constructs a spin foam model for this simplicial theory 
by imposing the constraints Vt %3 directly on the SU(2) BF amplitudes. The spin foam 
path integral for BF theory is obtained as in Section 0] The constraints are imposed by 
promoting the B l to operators x % (combinations of left / right invariant vector fields) acting 
on the discrete connection 17 . The model is defined as 

JV[ip] e^ii"" 

Zgr = I II d 9e 8(fF) E II A >f Tr IM ■ ■ ■ 9e)] . (46) 

c-.{j}^{f} feJ A 

^ .. ' 



JV[A] J V[B] e iJ Tr[BAF(A)] 

where Cl = J 1 A J 3 — |5* 3 J h A Jk and we have indicated the correspondence of the 
different terms with the continuum formulation. The preceding equation is rather formal; 



17 Notice that (for example) the right invariant vector field J l {U) — a 1 ^U B c d j dU A c has a well defined 
action at the level of equation (|32|l and acts as a B operator at the level of 129(1 since 

- i X \U) [e lTr [ BC/ l] \ v „ x = Tria'UB^^lu^ ~ B l ^ BU \ (45) 

where <r l are Pauli matrices. 
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for the rigorous implementation see j^Zj- Reisenberger uses locality so that constraints 
are implemented on a single 4-simplex amplitude. There is however a difficulty with the 
this procedure: the algebra of operators f2 u do not close so that imposing the constraints 
sharply becomes a too strong condition on the BF configurations 18 . In order to avoid this, 
Reisenberger defines a one-parameter family of models by inserting the operator 



2z 2 



(47) 



instead of the delta function above. In the limit z — > oo the constraints are sharply 
imposed. This introduces an extra parameter to the model. The properties of the kernel 
of Q have not been studied in detail. 



5.2 The Freidel-Krasnov prescription 

In reference [7T] Freidel and Krasnov define a general framework to construct spin foam 
models corresponding to theories whose action has the general form 

S(B, A) = J Tr [B A F(A)} + $(5), (48) 

where the first term is the BF action while $(-B) is a certain polynomial function of the B 
field. The formulation is constructed for compact internal groups. The definition is based 
on the formal equation 

I V[B}V[A] e ifn[BAF(A)]+HB) e ifH£) Z [J] _ , (49) 
where the generating functional Z[J) is defined as 

Z[J\ := / V[B]V[A] e </1>[BAF(A)]+lY[BAJ] ) (5Q) 



where J is an algebra valued 2-form field. They provide a rigorous definition of the gener- 
ating functional by introducing a discretization of M. in the same spirit of the other spin 
foam models discussed here. Their formulation can be used to describe various theories of 
interest such as BF theories with cosmological terms, Yang-Mills theories (in 2 dimensions) 
and Riemannian self-dual gravity. In the case of self dual gravity B and A are valued in 
su(2), while 



$(£) = / if, 



B l A B J - -5 lJ B k A B k 
3 



(51) 



according to equation (J43)) . The model obtained in this way is very similar to Reisenberger's 
one. There are however some technical differences. One of the more obvious one is that 
the non- commutative invariant vector fields J 1 representing B 1 are replaced here by the 
commutative functional derivatives 5/5 J' 1 . The explicit properties of these models have 
not been studied further. 



3 This difficulty also arises in the Barrett-Crane model as we shall see in Section 
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5.3 The Iwasaki model 



Iwasaki defines a spin foam model of self dual Riemannian gravity by a direct lattice 
discretization of the continuous Ashtekar formulation of general relativity. This model 
constitutes an example of the general prescription of Section 13.21 The action is 

S(e, A) = J dx' [2 + e° ljk e^ k u ] [2 d [x A^ + e°% m A l x A™} , (52) 

where A % a is an SU(2) connection. The fundamental observation of is that one can write 
the discrete action in a very compact form if we encode part of the degrees of freedom of 
the tetrad in an SU{2) group element. More precisely, if we take = e 1 e?J>ij = 1 we can 
define e M := e° + i<J^ G SU(2) where <jj are the Pauli matrices. In this parameterization 
of the 'angular' components of the tetrad and using a hypercubic lattice the discrete action 
becomes 

Sa = -P ^ Xa V» Tr H^U Xa ] , (53) 
veA 

where r M := (P 1 ^ 2 i P )~ 1 ey^gJ^l, U^ u is the holonomy around the /zz/-plaquette, e the lattice 
constant and (3 is a cutoff for r M used as a regulator (r M < /J 1 / 2 ^ -1 ). The lattice path 
integral is defined by using the Haar measure both for the connection and the 'spherical' 
part of the tetrad e's and the radial part dr^ := dr^r^. The key formula to obtain an 
expression involving spin foams is 

e ^^(2 i + l)^ Xi ([/). (54) 
J 

Iwasaki writes down an expression for the spin foam amplitudes in which the integration 
over the connection and the e's can be computed explicitly. Unfortunately, the integration 
over the radial variables r involves products of Bessel functions and its behavior is not 
analyzed in detail. In 3 dimensions the radial integration can be done and the corresponding 
amplitudes coincide with the results of Section (|4.3|) . 

5.4 The Barrett-Crane model 

The appealing feature of the previous models is the clear connection to loop quantum 
gravity, since they are defined directly using the self dual formulation of gravity (boundary 
states are S , f/(2)-spin networks). The drawback is the lack of closed simple expressions 
for the amplitudes which complicates their analysis. There is however a simple model 
that can be obtained as a systematic quantization of simplicial SO (4) Plebanski's action. 
This model was introduced by Barrett and Crane in [H] and further motivated by Baez in 
[Ti] . The basic idea behind the definition was that of the quantum tetrahedron introduced 

19 Iwasaki defines another model involving multiple cellular complexes to provide a simpler representation 
of wedge products in the continuum action. A more detail presentation of this model would require the 
introduction of various technicalities at this stage so we refer the reader to 72 . 
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by Barbieri in [7^ and generalized to 4d in ^5]. The beauty of the model resides in its 
remarkable simplicity. This has stimulated a great deal of explorations and produced many 
interesting results. We will review most of these in Section 

5.5 Markopoulou-Smolin causal spin networks 

Using the kinematical setting of LQG with the assumption of the existence of a micro-local 
(in the sense of Planck scale) causal structure Markopoulou and Smolin define a general 
class of (causal) spin foam models for gravity [1211101 (see also |77j). The elementary tran- 
sition amplitude A Sl ^ S[+1 from an initial spin network si to another spin network sj + i is 
defined by a set of simple combinatorial rules based on a definition of causal propagation of 
the information at nodes. The rules and amplitudes have to satisfy certain causal restric- 
tions (motivated by the standard concepts in classical Lorentzian physics). These rules 
generate surface-like excitations of the same kind we encounter in the more standard spin 
foam model but endow the foam with a notion of causality. Spin foams J-^^g, are labeled 
by the number of times these elementary transitions take place. Transition amplitudes are 



The models are not related to any continuum action. The only guiding principles are 
the restrictions imposed by causality, simplicity and the requirement of the existence of a 
non-trivial critical behavior that would reproduce general relativity at large scales. Some 
indirect evidence of a possible non trivial continuum limit has been obtained in some 
versions of the model in 1 + 1 dimensions. 

5.6 Gambini-Pullin model 

In reference [THJ Gambini and Pullin introduced a very simple model obtained by modifi- 
cation of the BF theory skein relations. As we argued in Section 0] skein relations defining 
the physical Hilbert space of BF theory are implied by the spin foam transition amplitudes. 
These relations reduce the large kinematical Hilbert space of BF theory (analogous to that 
of quantum gravity) to a physical Hilbert space corresponding to the quantization of a fi- 
nite number of degrees of freedom. Gambini and Pullin define a model by modifying these 
amplitudes so that some of the skein relations are now forbidden. This simple modification 
frees local excitations of a field theory. A remarkable feature is that the corresponding 
physical states are (in a certain sense) solutions to various regularizations of the scalar 
constraint for (Riemannian) LQG. The fact that physical states of BF theory solve the 
scalar constraint is well known [21], since roughly F(A) = implies EEF(A) = 0. The 
situation here is of a similar nature, and — as the authors argue — one should interpret this 
result as an indication that some 'degenerate' sector of quantum gravity might be repre- 
sented by this model. The definition of this spin foam model is not explicit since the theory 
is directly defined by the physical skein relations. 



defined as 




(55) 



N 1=0 
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5.7 Capovilla-Dell-Jacobson theory on the lattice 

The main technical difficulty that we gain in going from 3- dimensional general relativ- 
ity to the 4-dimensional one is that the integration over the e's becomes intricate. In 
the Capovilla-Dell-Jacobson |%T)| IHT] formulation of general relativity this 'integration' is 
partially performed at the continuum level. The action is 



where e ■ F A F := e F^F^. Integration over rj can be formally performed in the path 
integral and we obtain 



This last expression is of the form (jl7j) and can be easily discretized along the lines of 
Section 13.21 The final expression (after integrating over the lattice connection) involves 
a sum over spin configurations with no implicit integrations. One serious problem of this 
formulation is that it corresponds to a sector of gravity where the Weyl tensor satisfy certain 
algebraic requirements. In particular flat geometries are not contained in this sector. 

6 Some conceptual issues 
6.1 Anomalies and gauge fixing 

As we mentioned before and illustrated with the example of Section the spin foam path 
integral is meant to provide a definition of the physical Hilbert space. Spin foam transition 
amplitudes are not interpreted as defining propagation in time but rather as defining the 
physical scalar product. This interpretation of spin foam models is the only one consistent 
with general covariance. However, in the path integral formulation, this property relies on 
the gauge invariance of the path integral measure. If the measure meets this property we 
say it is anomaly free. It is well known that in addition to the invariance of the measure, 
one must provide appropriate gauge fixing conditions for the amplitudes to be well defined. 
In this section we analyze these issues in the context of the spin foam approach. 

Since we are interested in gravity in the first order formalism, in addition to diffeomor- 
phism invariance one has to deal with the gauge transformations in the internal space. Let 
us first describe the situation for the latter. If this gauge group is compact then anomaly 
free measures are defined using appropriate variables and invariant measures. In this case 
gauge fixing is not necessary for the amplitudes to be well defined. Examples where this 
happens are: the models of Riemannian gravity considered in this paper (the internal gauge 
group being £0(4) or SU(2)), and standard lattice gauge theory. In these cases, one rep- 
resents the connection in terms of group elements (holonomies) and uses the (normalized) 
Haar measure in the integration. In the Lorentzian sector (internal gauge group SL(2, C)) 




(56) 




X 



(57) 
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the internal gauge orbits have infinite volume and the lattice path integral would diverge 
without an appropriate gauge fixing condition. These conditions generally exist in spin 
foam models and we will study an example in Section El (for a general treatment see [82]). 

The remaining gauge freedom is diffeomorphism invariance. It is generally assumed 
that spin foam encode diff-invariant information about 4-geometries. Thus, no gauge fix- 
ing would be necessary as one would be already summing over physical configurations. 
To illustrate this perspective we concentrate on the case of a model defined on a fixed 
discretization A as described in Section l3~2l 

Let us start by considering the spin network states, at <9A: boundary of A, for which 
we want to define the transition amplitudes. According to what we have learned from 
the canonical approach, 3-diffeomorphism invariance is implemented by considering (dif- 
feomorphism) equivalence classes of spin-network states. In the context of spin foams, the 
underlying discretization A restricts the graphs on the boundary to be contained on the 
dual 1-skeleton of the boundary complex dA. These states are regarded as representative 
elements of the corresponding 3-diffeomorphism equivalence class. The discretization can 
be interpreted, in this way, as a gauge fixing of 3-diffeomorphisms on the boundary. This 
gauge fixing is partial in the sense that, generically, there will remain a discrete symmetry 
remnant given by the discrete symmetries of the spin network. This remaining symmetry 
has to be factored out when computing transition amplitudes (in fact this also plays a role 
in the definition of the kinematical Hilbert space of Section EJ). 

The standard view point (consistent with LQG and quantum geometry) is that this 
should naturally generalize to 4-difleomorphisms for spin foams. The underlying 2-complex 
J7a on which spin foams are defined represents a partial gauge fixing for the configurations 
(spin foams) entering in the path integral. The remaining symmetry, to be factored out in 
the computation of transition amplitudes, corresponds simply to the finite group of discrete 
symmetries of the corresponding spin foams 20 . This factorization is well defined since the 
number of equivalent spin foams can be characterized in a fully combinatorial manner, 
and is finite for any spin foam defined on a finite discretization. In addition, a spin foam 
model is anomaly free if the amplitudes are invariant under this discrete symmetry. This 
requirement is met by all the spin foam models we considered in this article. 

We illustrate the intuitive idea of the previous paragraph in Figure El On the diagram 
(a) a continuum configuration is represented by a discrete (spin foam) configuration on 
the lattice. On (b) two configurations are shown. In the background dependent context 
(e.g., lattice gauge theory) these two configurations would be physically inequivalent as 
the lattice carries metric information (length of the edges). In the context of spin foam 
models there is no geometric information encoded in the discretization and in an anomaly 
free spin foam model the two configurations should be regarded as equivalent. 

The discretization of the manifold A is seen as a regulator introduced to define the 
spin foam model. Even when the regulator (or the discretization dependence) eventually 
has to be removed (see next subsection), the theory remains discrete at the fundamental 

20 Baez |14j points out this equivalence relation between spin foams as a necessary condition for the 
definition of the category of spin foams. 
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Figure 6: Diffeomorphisms are replaced by discrete symmetries in spin foams. In this 
simplified picture there are two possible colorings of or 1, the latter represented by the 
darkening of the corresponding lattice element. The two configurations on (b) are regarded 
as physically equivalent in the background independent context. 



level. The smooth manifold diffeomorphism invariant description is regarded in this context 
as a derived concept in the (low energy) continuum limit. Fundamental excitations are 
intrinsically discrete. From this viewpoint, the precise meaning of the gauge symmetries 
of such a theory would have to be formulated directly at the discrete level. We have seen 
that this can be achieved in the case of 3-dimensional gravity (recall Section 14.3.1)) . From 
this perspective, the discrete symmetries of colored 2-complexes (spin foams) represent the 
fundamental 'gauge' freedom. This symmetry would manifest itself as diffeomorphisms 
only in the continuum limit. 

The viewpoint stated above is consistent with recent results obtained by Gambini and 
Pullin |83| IH4"] . They study the canonical formulation of theories defined on a lattice from 
the onset. This provides a way to analyze the meaning of gauge symmetries directly a 
la Dirac. Their results indicate that diff-invariance is indeed broken by the discretization 
in the sense that there is no infinitesimal generator of diffeomorphism. This is consistent 
with our covariant picture of discrete symmetries above. In their formulation the canonical 
equations of motion fix the value of what were Lagrange multipliers in the continuum (e.g. 
lapse and shift). This is interpreted as a breaking of diffeomorphism invariance; however, 
the solutions of the multiplier equations are highly non unique. The ambiguity in selecting 
a particular solution corresponds to the remnant diffeomorphism invariance of the discrete 
theory. 

However, the notion of anomaly freeness stated at the beginning of this section should 
be strengthened. In fact according to our tentative definition, an anomaly free measure 
can be multiplied by any gauge invariant function and yield a new anomaly free measure. 
This kind of ambiguity is not what we want; however, it is in fact present in most of the 
spin foam models defined so far. In standard QFT theory, the formal (phase space) path 
integral measure in the continuum has a unique meaning (up to a constant normaliza- 
tion) emerging from the canonical formulation. Provided an appropriate gauge fixing, the 
corresponding Dirac bracket determines the formal measure on the gauge fixed constraint 
surface. As we will see later, there is a certain degree of ambiguity in the definition of 
various spin foam models in four dimensions. This ambiguity concerns the evaluation of 
lower dimensional simplexes and is directly related to the definition of measure in the spin 
foam sum. One would expect that a strengthened definition of anomaly freeness should 
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resolve these ambiguities. This possibility is studied in 

Finally let us briefly recall the situation in 3- dimensional gravity. In three dimensions 
the discrete action is invariant under transformations that are in correspondence with the 
continuum gauge freedoms (f2*3*|) and (|24jl . As we mentioned in Section [4.3.11 spin foam 
amplitudes diverge due to (J2IJ). This means that spin foams do not fix (f2"4"j) up to a 'finite 
volume' discrete symmetry, which seems to be in conflict with the argument stated above as 
diffeomorphism can be obtained combining the transformations ()23|) and (|24p. However, 
the topological gauge symmetry ()24j) involves more than just diffeomorphisms. Indeed, 
only on shell one can express diffeomorphisms as a combination of (|2*3*|) and ()24jl . This 
representation of diffeomorphisms turns out to be field dependent (recall Footnote ITTj) . At 
the canonical level we have seen that the topological symmetry (|24|) acts in away that is 
totally different from diffeomorphisms. In particular spin network states defined on graphs 
which which differ in the number of edges and vertices can be physically equivalent (recall 
Section T4.3.3p . Thus 3-dimensional gravity is a degenerate example and extrapolation of 
its gauge properties to four dimensions seems misleading. 

In Section 17.11 we will study the spin foam quantization of BF theory — a topological 
theory which corresponds to a generalization of 3-dimensional gravity to four dimensions. 
The spin foam amplitudes are also divergent in this case due to the analog of the gauge 
symmetry (j24j). Some of the spin foam models in Section El are defined from BF theory 
by implementing constraints that reduce the topological theory to general relativity. The 
implementation of the constraints breaks the topological symmetry (J2*4"|) and the resulting 
model is no longer topological invariant. At the continuum level it is clear that the remnant 
gauge symmetry is diffeomorphism invariance. Whether in the the resulting spin foam 
model for gravity remnant is larger than the discrete spin foam symmetries advocated 
above has not been studied in detail. 

The action of diffeomorphism is far from understood at a rigorous level in the context 
of spin foam models. Here we have presented an account of some of the ideas under 
consideration, and we tried to point out the relevance of an issue that certainly deserves 
detailed investigation. 

6.2 Discretization dependence 

The spin foam models we have introduced so far are defined on a fixed cellular decom- 
position of M.. This is to be interpreted as an intermediate step toward the definition 
of the theory. The discretization reduces the infinite dimensional functional integral to a 
multiple integration over a finite number of variables. This cutoff is reflected by the fact 
that only a restrictive set of spin foams (spin network histories) is allowed in the path 
integral: those that can be obtained by all possible coloring of the underlying 2-complex. 
In addition it restricts the number of possible 3-geometry states (spin network states) on 
the boundary by fixing a finite underlying boundary graph. This represents a truncation 
in the allowed fluctuations and the set of states of the theory that can be interpreted 
as a regulator. However, the nature of this regulator is fundamentally different from the 
standard concept in the background independent framework: since geometry is encoded in 



31 



the coloring (that can take any spin values) the configurations involve fluctuations all the 
way to Plank scale 21 . This scenario is different in lattice gauge theories where the lattice 
introduces an effective UV cutoff given by the lattice spacing. Transition amplitudes are 
however discretization dependent now. A consistent definition of the path integral using 
spin foams should include a prescription to eliminate this discretization dependence. 

A special case is that of topological theories such as gravity in 3 dimensions. In this 
case, one can define the sum over spin foams with the aid of a fixed cellular decomposition 
A of the manifold. Since the theory has no local excitations (no gravitons), the result is 
independent of the chosen cellular decomposition. A single discretization suffices to capture 
the degrees of freedom of the topological theory. 

In lattice gauge theory the solution to the problem is implemented through the so-called 
continuum limit. In this case the existence of a background geometry is crucial, since it 
allows one to define the limit when the lattice constant (length of links) goes to zero. In 
addition the possibility of working in the Euclidean regime allows the implementation of 
statistical mechanical methods. 

None of these structures are available in the background independent context. The 
lattice (triangulation) contains only topological information and there is no geometrical 
meaning associated to its components. As we mentioned above this has the novel con- 
sequence that the truncation can not be regarded as an UV cutoff as in the background 
dependent context. This in turn represents a conceptual obstacle to the implementation 
of standard techniques. Moreover, no Euclidean formulation seems meaningful in a back- 
ground independent scenario. New means to eliminate the truncation introduced by the 
lattice have to be developed. 

This is a major issue where concrete results have not been obtained so far beyond the 
topological case. Here we explain the two main approaches to recover general covariance 
corresponding to the realization of the notion of 'summing over discretizations' of p>f] . 

• Refinement of the discretization: 

According to this idea topology is fixed by the simplicial decomposition. The trun- 
cation in the number of degrees of freedom should be removed by considering tri- 
angulations of increasing number of simplexes for that fixed topology. The flow in 
the space of possible triangulations is controlled by the Pachner moves. The formal 
idea is to take a limit in which the number of four simplexes goes to infinity together 
with the number of tetrahedra on the boundary. Given a 2-complex Ji which is a 
refinement of a 2-complex J\ then the set of all possible spin foams defined on J\ 
is naturally contained in those defined on J 2 (taking into account the equivalence 
relations for spin foams mentioned in the previous section). The refinement process 
should also enlarge the space of possible 3-geometry states (spin networks) on the 
boundary recovering the full kinematical sector in the limit of infinite refinements. 
An example where this procedure is well defined is Zapata's treatment of the Turaev- 
Viro model [70] . The key point in this case is that amplitudes are independent of the 

21 Changing the label of a face from j to j + 1 amounts to changing an area eigenvalue by an amount of 
the order of Planck length squared according to Q . 
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discretization (due to the topological character of the theory) so that the refinement 
limit is trivial. In the general case there is a great deal of ambiguity involved in the 
definition of refinement 22 . The hope is that the nature of the transition amplitudes 
would be such that these ambiguities will not affect the final result. The Turaev-Viro 
model is an example where this prescription works. 

If the refinement limit is well defined one would expect that working with a 'suffi- 
ciently refined' but fixed discretization would serve as an approximation that could 
be used to extract physical information with some quantifiable precision. 

• Spin foams as Feynman diagrams: 

This idea has been motivated by the generalized matrix models of Boulatov and 
Ooguri [SHIEZI- The fundamental observation is that spin foams admit a dual for- 
mulation in terms of a field theory over a group manifold ^IJ EJ EI- The duality 
holds in the sense that spin foam amplitudes correspond to Feynman diagram ampli- 
tudes of the GFT. The perturbative Feynman expansion of the GFT (expansion in a 
fiducial coupling constant A) provides a definition of sum over discretizations which 
is fully combinatorial and hence independent of any manifold structure 23 . The latter 
is most appealing feature of this approach. 

However, the convergence issues clearly become more involved. The perturbative 
series are generically divergent. This is not necessarily a definite obstruction as 
divergent series can often be given an asymptotic meaning and provide physical 
information. Moreover, there are standard techniques that can allow to 're-sum' a 
divergent series in order to obtain non perturbative information. Recently, Freidel 
and Louapre jH5] have shown that this is indeed possible for certain GFT's in three 
dimensions. Other possibilities have been proposed in [To] . 

It is not clear how the notion of diffeomorphism would be addressed in this framework. 
Diffeomorphism equivalent configurations (in the discrete sense described above) ap- 

22 It is not difficult to define the refinement in the case of a hypercubic lattice. In the case of a simplicial 
complex a tentative definition can be attempted using Pachner moves. To illustrate we can concentrate on 
the simple 2-dimensional case. Given an initial triangulation of a surface with boundary Ai define Aj by 
implementing a 1 — 3 Pachner move to each triangle in Ai . A^ is a homogeneous refinement of Ai ; however, 
the boundary triangulation remain unchanged so that they will support the same space of boundary data 
or spin networks. As mentioned above, in the refinement process one also wants to refine the boundary 
triangulation so that the corresponding dual graph will get refined and the space of possible boundary 
data will become bigger (it should involve all of Cyl in the limit). In order to achieve this we define A 2 by 
erasing from Aj triangles sharing a 1-simplex with the boundary. This amounts to carrying out a 1 — 2 
Pachner move on the boundary. This completes the refinement Ai — > A 2 . This refinement procedure 
seems however not fully satisfactory as the refinement keeps the memory of the initial triangulation. This 
can be easily visualized in the dual 2-complex. An improvement of this prescription could involve the 
random implementation of 2 — 2 Pachner moves between refinements. 

23 This is more than a 'sum over topologies' as many of the 2-complex appearing in the perturbative 
expansion cannot be associated to any manifold [BB| . 
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pear at all orders in the perturbation series 24 . From this perspective (and leaving 
aside the issue of convergence) the sum of different order amplitudes corresponding 
to equivalent spin foams should be interpreted as the definition of the physical ampli- 
tude. The discussion of the previous section does not apply in the GFT formulation, 
i.e., there is no need for gauge fixing. 

The GFT formulation resolves by definition the two fundamental conceptual prob- 
lems of the spin foam approach: diffeomorphism gauge symmetry and discretization 
dependence. The difficulties are shifted to the question of the physical role of A and 
the convergence of the corresponding perturbative series. The GFT formulation has 
also been very useful in the definition of a simple normalization of the Barret-Crane 
model and has simplified its generalizations to the Lorentzian sector. We will study 
this formulation in detail in Section |HJ 



6.3 Physical scalar product revisited 

In this subsection we study the properties of the generalized projection operator P intro- 
duced before. The generalized projection operator is a linear map from the kinematical 
Hilbert space into the physical Hilbert space (states annihilated by the scalar constraint). 
These states are not contained in the kinematical Hilbert space but are rather elements of 
the dual space Cyl*. For this reason the operator P 2 is not well defined. We have briefly 
mentioned the construction of the generalized projection operator in the particular context 
of Rovelli's model of Section (jfl.lj) . After the above discussion of gauge and discretization 
dependence we want to revisit the construction of the P operator. We consider the case in 
which discretization independence is obtained by a refinement procedure. 

The number of 4-simplexes N in the triangulation plays the role of cutoff. The matrix 
elements of the P operator is expected to be defined by the refinement limit of fixed 
discretization transition amplitudes ((sPn, s')), namely 

(*, s') phys = (sP, s') = lim (sP N , s') . (58) 

If we have a well defined operator P (formally corresponding to 5[S]) the reconstruction 
of the physical Hilbert space goes along the lines of the GNS construction j^I] where the 
algebra is given by Cyl and the state ufif) for \P e Cyl is defined by P as 

«<»> = 

V ; (OP, 0) ' 

In the spin foam language this corresponds to the transition amplitude from the state \I/ 
to the 'vacuum'. The representation of this algebra implied by the GNS construction has 
been used in [92J to define observables in the theory. 

24 The GFT formulation is clearly non trivial already in the case of topological theories. There has been 
attempts to make sense of the GFT formulation dual to BF theories in lower dimensions 
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In Section 13.11 we mentioned another way to obtain observables according to the defi- 
nition 

(sOphys, s } = (sPO kin P, s ), (59) 

where O kin is a kinematical observable. This construction is expected to be well defined 
only for some suitable kinematical operators O kin (e.g., the previous equation is clearly 
divergent for the kinematical identity). There are kinematical operators which can be used 
as regulators 25 . They corresponds to Rovelli's 'sufficiently localized in time' operators [53] . 
In [HSj it is argued that in posing a physical question one is led to conditioning the path 
integral and in this way improving the convergence properties. Conditioning is represented 
here by the regulator operators (see previous footnote). 

25 Consider the following simple example. Take the torus S 1 x S 1 represented by [0,27r] x [0,27r] € R 2 
with periodic boundary conditions. An orthonormal basis of states is given by the ('spin network') wave 
functions (xy,nm) = J- e mx e lm v for n,m e Z. Take the scalar constraint to be given by X so that the 
generalized projection P simply becomes 
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The RHS is no longer in Jikin since is not normalizable: it is meaningful as a distributional state. The 
physical scalar product is defined as in (|58(l . namely 

(nm,n'm') phys = (nmP,n'm') = ^ (am, n'm') = S mm >. 

OL 

Therefore we have that \nm) = \km) Vfe. We can define physical observables by means of the formula 

(m, O p h m') — (nm, POP n'm') , 

for suitable kinematic observables O. If for example we take = 1 the previous equation would involve 
P 2 which is not defined. However, take the following family of kinematical observables 

\O k ,nm) = 5 n> k | fern) , 

then (rn,Oph m ') — S m .m', i.e., O p h is just the physical identity operator for any value of k. The whole 
family of kinematic operators define a single physical operator. The key of the convergence of the previous 
definition is in the fact that Ok is a projection operator into an eigenstate of the operator P x canonically 
conjugate to the constraint X. In particular any operator that decays sufficiently fast along going to ±oo in 
the spectrum of p x will do. This illustrates what kind of operators will be suitable for the above definition 
to work: they have to be sufficiently localized along the gauge orbit, hence 'sufficiently localized in time' 

— (fc — n)^ 

in the case of the scalar constraint. For example, one may also use the Gaussian Ok \nm) = e i \nm) 
which also projects down to a multiple of the physical identity. Moreover, availability of such a kinematical 
operator can serve as a regulator of non suitable operators, for instance 



= (nm,Pd k dO k P n'm'), 



7 Oph 

will converge even for = 1. The operator Ok can be regarded as a quantum gauge fixing: it selects a 
point along the gauge orbit generated by the constraint. 
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6.4 Contact with the low energy world 

A basic test to any theory of quantum gravity is the existence of a well defined classical 
limit corresponding to general relativity. In the case of the spin foam approach this should 
be accomplished simultaneously with a limiting procedure that bridges the fundamental 
discrete theory with the smooth description of classical physics. That operation is some- 
times referred to as the continuum limit but it should not be confused with the issues 
analyzed in the previous subsection. 

There is debate on how one would actually setup the definition of the 'low energy 
limit' in the background independent context. The general strategy — in fact motivated by 
our experience in background physics — is to setup a renormalization scheme (a la Wilson) 
where microscopic degrees of freedom are summed out to obtain a coarse grained effective 
description. The previous strategy is to be regarded as heuristic until a clear-cut definition 
is found which encompasses all the issues discussed in this section. It should be emphasized 
that even when various spin foam models for quantum gravity have been defined, none of 
them has been show to reproduce gravity at 'low energies'. This is a major open question 
that deserves all our efforts. 

Interesting ideas on how to define the renormalization program for spin foams have 
been proposed by Markopoulou jHH EE] • They make use of novel mathematical techniques 
shown to be useful in dealing with Feynman diagrams in standard QFT Wf\ . In this 
approach Pachner moves on the underlying 2-complex are regarded as radiative corrections 
in standard QFT 26 . In order to understand the 'scaling' properties of the theory one has to 
study the behavior of amplitudes under the action of these elementary moves. Amplitudes 
for this elementary moves in the Barrett-Crane model have been computed in jl2j. A full 
implementation of Markopoulou's proposal to a concrete model remains to be studied. 

Recently Oeckl [2H] has studied the issue of renormalization in the context of spin foam 
models containing a coupling parameter. These models include generalized covariant gauge 
theories P^lllUUj . the Reisenberger model, and the so called interpolating model (defined by 
Oeckl). The latter is given by a one-parameter family of models that interpolate between 
the trivial BF topological model and the Barrett-Crane model according to the value of 
a 'coupling constant'. Qualitative aspects of the renormalization groupoid flow of the 
couplings are studied in the various models. 



26 this interpretation becomes completely precise in the GFT formulation of spin foams presented in 
Section El 
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Part II 



The Barrett-Crane model 

7 SO (4) Plebanski's action and the Barrett-Crane model 

The Barrett-Crane model is one of the most extensively studied spin foam models for 
quantum gravity. In this section we concentrate on the definition of the model in the 
Riemannian sector and sketch its generalization to the Lorentzian sector. 

The Barrett-Crane model can be viewed as a spin foam quantization of SO (4) Pleban- 
ski's formulation of general relativity. The general idea is analogous to the one implemented 
in the model of Section IHTT1 We introduce the model from this perspective in the Rieman- 
nian sector so we need to start with the study of Spin(A) BF theory. 

7.1 Quantum Spin (4) BF theory 

Classical (Spin(A)) BF theory is defined by the action 

S[B, A] = J Ti[B A F(A)} , (60) 

M 

where B^ is a Spin(4) Lie-algebra valued 2-form, Aj/ is a connection on a Spin(4) principal 
bundle over M.. The theory has no local excitations. Its properties are very much analogous 
to the case of 3-dimensional gravity studied in Section 0J 

A discretization A of M. can be introduced along the same lines presented in Section 
@] For simplicity we concentrate in the case when A is a triangulation. The field B is 
associated with Lie algebra elements B t assigned to triangles t £ A. In four dimensions 
triangles t £ A are dual to faces / £ J7a- This one-to-one correspondence allows us to 
denote the discrete B by either a face (£>/) or a triangle (B t ) subindex respectively. Bf 
can be interpreted as the 'smearing' of the continuous 2-form B on triangles in A. The 
connection A is discretized by the assignment of group elements g e £ Spin(4) to edges 
e £ c7a- The path integral becomes 

Z(A) = f J] dg e J] dBf e iB f u f = f \{ dg e J] 5(g ei ■ ■ ■ g e J, (61) 

where the first equality is the analog of (|2T?|) while the second is the result of the B inte- 
gration |101j . Using Peter- Weyl's theorem as in (p?2"j) one obtains 

Z(A) = l II d 9e II A Pf Tr IPfti ■ ■ ■ 9e)] , (62) 

c--{p}->{f} J eej A feJ A 

where p are Spin(A) irreducible representations. Integration over the connection can be 
performed as in the case of 3-dimensional gravity. In a triangulation A, the edges e £ J7a 
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bound precisely four different faces; therefore, the g e 's in (R)2*)l appear in four different 
traces. The relevant formula is 

Pi v := I dg Pi{g)®P2{g)®---®PA{g) = Y, C n P ^ C Zn--^ ( 63 ) 

where Pf nv is the projector onto Inv[pi ® p 2 ® • • • <8> p n ] and on the RHS we have expressed 
the projector in terms of normalized intertwiners as in (J33|) . Finally, integration over the 
connection yields 
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zbf(a)= £ yi n^n 

C/:{/}-*/3/C e :{e}->{te} /£Ja v£j A 




where the pentagonal diagram representing the vertex amplitude denotes the trace of the 
product of five normalized 4-intertwiners C L . As in the model for 3-dimensional 

1 P1P2P3P4 

gravity, the vertex amplitude corresponds to the flat evaluation of the spin network state 
defined by the pentagonal diagram in (|64jh a 15j-symbol. Vertices v G J7a are in one-to-one 
correspondence to 4-simplexes in the triangulation A. In addition we have C e : {e} — > {i e } 
representing the assignment of intertwiners to edges coming from (I63|) . 

The state sum is generically divergent (due to the gauge freedom analogous to 
f|24)l ). A regularized version defined in terms of SU q (2) x SU q {2) was introduced by Crane 
and Yetter 1102j . As in three dimensions is topologically invariant and the spin 
foam path integral is discretization independent. 



7.2 Classical 50(4) Plebanski action 

Plebanski's Riemannian action can be thought of as Spin(4) BF theory plus constraints 
on the S-field. It depends on an so(4) connection A, a Lie-algebra-valued 2-form B and 
Lagrange multiplier fields A and p. Writing explicitly the Lie-algebra indexes, the action 
is given by 

S[B, A, A,/i] = J [B IJ A Fjj{A) + Xukl B IJ A B KL + pe IJKL \ IJKL ] , (65) 

where \i is a 4-form and Xukl = —Xjikl = —Xulk = Xklij is a tensor in the internal 
space. Variation with respect to // imposes the constraint e 13KL \ijKL = on Xukl- The 
Lagrange multiplier tensor Xukl has then 20 independent components. Variation with 
respect to A imposes 20 algebraic equations on the 36 components of B. Solving for \x they 
are 

e »vpa B IJ B KL = e e IJKL (gg) 



38 



which is equivalent to 



£ijklB*IB x = e e^ upa , (67) 



pv pa 

for e 7^ where e = ^opqrB?^ BQ^e llvp<J jlU3j . The solutions to these equations are 

B = ±*(eAe), and B = ±e A e, (68) 

in terms of the 16 remaining degrees of freedom of the tetrad field e a . If one substitutes the 
first solution into the original action one obtains Palatini's formulation of general relativity 

S[ e , A] = J Tr [e A e A *F{A)\ . (69) 

This property is the key to the definition of the spin foam model for gravity of the next 
section. 



7.3 Discretized Plebanski's constraints 

In Section 17.11 we have derived the spin foam model for Spin (4) BF theory. We have 
seen that the sum over ^-configurations is encoded in the sum over unitary irreducible 
representations of Spin(A). Can we restrict the spin foam configurations to those that 
satisfy Plebanski's constraints? To answer this question one first needs to translate the 
constraints of Plebanski's formulation into the simplicial framework. There are two ways 
in which we can approach the problem. One is by looking at the constraints in the form 
(J55J) the other is by using (JH7J). 

The first version of the constraints is of the type of of Section 15.11 This possibility 
has not been explored in the literature. A simple representation of the constraints is 
obtained if we choose to discretize (|67|). for which there are no free algebra indexes. 

As explained in Section I7~T1 in the discretization of the B field we replace the space-time 
indexes of the 2-form by a triangle: 

iw — > t = triangle (70) 

Using this it is easy to translate (|67|) into the simplicial framework. The RHS of (J67|) 
vanishes when there are repeated space-time indexes. At the discrete level, this means 
that whenever the triangles t and t' belong to the same tetrahedron T, t,t' £ T, (i.e., 
either t — t' or t and t' share an edge) we have 

e IJKL Bl J B^ L = 0, (71) 

The remaining constraints correspond to the situation when t and t' do not lie on the same 
tetrahedron (no 'repeated space-time indexes'). We can write these constraints as follows. 
If we label a vertex in a given 4-simplex, ■ ■ • 4 we can label the six triangles Ojk containing 
the vertex simply as jk (j, k — 1 • • - 4). Only these are needed to express the independent 
constraints 

cijklB^B^ 1 = euKiB^B^ 1 = euKLB^B^ 1 oc e. (72) 
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Now we want to implement the constraints. The strategy is to impose them directly 
on the BF spin foam sum, (|62jh after the 5-integration has been performed. Formally one 
would associate the discrete Bf to the differential operator —id/ dUf in (|ET|). More precisely, 
the observation is that the Spin(A) left invariant vector field —iX IJ (U) := £/ M U X IJ u a 9 ^ 
acts as a quantum B IJ on (JoTj) since 

-iX"(U) (e^) |^ = U\X IJ ^^L-e^l^ = 

= Ti[UX IJ B]e iT ^ BU %^ ~ B /J e ffl ™, (73) 

where X 7J are elements of an orthonormal basis in the 50(4) Lie- algebra. The evaluation 
at U — 1 is motivated by the fact that configurations in the BF partition function (jHTjl 
have support on flat connections. This approximation is made in order to motivate our 
definition but it plays no role in the implementation of the constraints. 

The constraints (JoTj) are quadratic in the B's. We have to worry about cross terms; the 
nontrivial case corresponds to: 

e IJKL X IJ {U)X KL {U) (e^l) |^ 

= - eiJKL (Tr[X IJ UB]Ti[X KL UB]e m ^ + iTi[X IJ X KL U B}^^) 
~e IJKL B IJ B KL ^ BU \ (74) 

where the second term on the second line can be dropped using that €uklX ij X oc 1 
(one of the two 50(4) Casimir operators) and U ~ 1. Therefore, we define the Bf field 
associated to a face at the level of equation (J32|) as the appropriate left invariant vector 
field —iX IJ (Uf) acting on the corresponding discrete holonomy Uf, namely 

B 1 / - -iX IJ (U f ). (75) 

Gauge invariance of the BF partition function implies that for every tetrahedron 

where t G T denotes the triangles in the corresponding tetrahedron. 

In order to implement (fTTj) we concentrate on a single 4-simplex amplitude, i.e., the 
fundamental atom of Figure 0] Once the constrained 4-simplex amplitude is constructed 
any spin foam amplitude can be obtained by gluing atoms together along faces by inte- 
gration over common boundary data as in Section 13.21 The BF 4-simplex wave function 
is obtained using (|61|) on the dual 2-complex with boundary defined by the intersection of 
the dual of a single 4-simplex with a 3-sphere, see Figure |3J 

The amplitude of the fundamental atom is a (cylindrical) function depending on the 
boundary values of the connection on the boundary graph 75. We denote as hij G 5pin(4) 
(i ^ j, i,j — 1 • • • 5 and = h~^) the corresponding 10 boundary variables (associated 
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to thin boundary edges in Figure 0} 27 and gi G Spin(A) (i = 1, • • • , 5) the internal connec- 
tion (corresponding to the thick edges in Figure EJ). According to (|61|1 the 4-simplex BF 
amplitude ASIMsFihij) is given by 

AsiM BF (h i3 ) = I n^n^^)' ( 7? ) 

where Uij = gih^Qj is the holonomy around the triangular face Oij. With the definition of 
the B fields given in (J75J) the constrained amplitude, 4S I M const (hij) , formally becomes 

ASIMconstihij) = / II d & 5 [Constraints(-zA'(Ui i ))] JJ 8^^). (78) 

i i<j 

It is easy to verify, using an equation analogous to (fTHj) and the invariance of €ukl, that 
one can define the B's by simply acting with the left invariant vector fields on the boundary 
connection hy. Therefore, the previous equation is equivalent to 

ASIM const (hij) = 5 [Constraints(-zAf(/ij i ))] / &(9ihijgj), (79) 

i i<j 

where we have taken the delta function out of the integral. The quantity on which the 
formal delta distribution acts is simply 4SIM BF (hij). The amplitude 4SIM BF (hij) can be 
expressed in a more convenient way if we expand the delta functions in modes as in (|62[) 
and then integrate over the internal connection g^. The integration is analogous to the one 
in (jSHJl) for example integration over g\ yields 




where on the RHS we have chosen a particular basis of 4-intertwiners to span the pro- 
jector Pf nv into Inv [p i2 ■ • - pis]. The 4-intertwiners are explicitly given as contractions of 
normalized intertwiners using the graphical notation introduced in Section 0J The circles 

27 Strictly speaking, the boundary connections hij are denned as the product h'^h^ where h! and h are 
associated to half paths as follows: take the edge ij for simplicity and assume it is oriented from i to j. 
Then h'^ is the discrete holonomy from i to some point in the center of the path and is the holonomy 
from that center point to j. This splitting of variables is necessary when matching different atoms to 
reconstruct the simplicial amplitude |57| . 
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represent the corresponding p-representation matrices evaluated on the the corresponding 
boundary connection h. The 4-simplex amplitude becomes 



la lJ 




In the previous equation 4-valent nodes denote normalized 4-intertwiners and the tree 
decomposition is left implicit (the factors A t in (jHOj) have been absorbed into the notation). 
The term on the left is a 15j-symbol as in ()64|) while the term on the right is the trace of 
five 4-intertwiners with the respective boundary connection insertions. 

The 4-simplex amplitude for the constrained spin foam model is then defined as the 
restriction of 4SIMBF(hij) imposed by the quantum version of the constraints (|7T|) defined 
by the following set of differential equations 

Cij, lk 4SIM const (h i3 ) = e IJKL X IJ (h l3 )X KL (h lk ) ASIM const (h l3 ) = V j, k, (82) 

and where the index i — 1, • • • , 5 is held fixed and Cy, it denotes the corresponding con- 
straint operator. 

There are six independent constraints (|H2*j) for each value of i — 1, • • • , 5. If we consider 
all the equations for the 4-simplex amplitude then some of them are redundant due to 
(I76|) . The total number of independent conditions is 20, in agreement with the number of 
classical constraints ([67)1 . For a given i in (|82j) (i.e., a given tetrahedron) and for j = k the 
equation becomes 

euK L X IJ (h ii )X KL {h i j) ASIM const {h l3 ) = 

= 5 ik [J^h^J^h^) - Xihl^ih^)] \SIM con MrK 3 ) = 

= + !) - fM + !)] 45/M const (^) = 0, (83) 

where we have used that: Spin{A) = SU{2) x SU{2) so that for h G Spin^i), h e ,ti r e 
SU(2) denote its right and left components, irreducible representation can be expressed as 
p = j e <8> j r for j e ,j r G Irrep[S77(2)], and the left invariant vector field 

X IJ (h) = P+ IJ J i (h t ) + P- Ij j\h r ), (84) 

for ejj KL P ±IJ = ±p ±KL and J l 's being left invariant vector fields on the corresponding 
left and right SU{2) copies of Spin{4). The previous constraints are solved by requiring 
the corresponding representation to be simple, i.e., 

pij = jij ® jij or Pij = Jij ® hi- (85) 
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This ambiguity is analogous to the classical one in (J55|) . We take p^- = jij ® j*j in corre- 
spondence to the choice *(e A e) that produces the gravity sector 28 . This solves 10 of the 
20 equations. The next non-trivial conditions imposed by (}82|) is when j ^ k. In this case 
we have 

2e IJKL X IJ (h ij )X KL (h lk ) ASIM const (h l3 ) 

= e IJKL (* /J (/^) + X IJ {h lk )) (X KL (h l3 ) + X KL (h lk )) ASIM const {h l3 ) 
= + 1) - L r (i r + 1)] ASIM const (h^) 

= 0, (86) 

where in the second line we used the fact that we have already solved (JBEJ)- In the third 
line we have used the gauge invariance (or the analog of (|7fi|) for the 3-valent node in the 
tree decomposition that pairs the representation p^ with the p ik ) in order to express the 
sum of invariant vector fields as the invariant vector filed acting on the virtual link labeled 
by l. This choice of tree decomposition in the case ij = 12 and ik = 13 is the one used in 
equation (JHUj) . The solution is clearly l — l ® l*. 

What happens now to any of the two remaining conditions (only one is independent), 
for example, C f y/ j w for k ^ k', j ^ j' and j' ^ k'l At first sight it looks like this equa- 
tions cannot be (generically) satisfied because an intertwiner that has simple i in one tree 
decomposition does not have only simple l' components in a different tree decomposition 
as a consequence of the recoupling identity (|42|l . There is however a linear combination of 
intertwiners found by Barrett and Crane in [Oj which is simple in any tree decomposition, 
namely 

\**o)= E ( 87 ) 

simple i 

C u ... p4 is a normalized 4-intertwiner and the summation is over simple i (i.e. i = i ® i*) 
and the pi are also simple (pi = ji ® j* for i = 1..4). This is clearly a solution to all 
the constraints and has been shown to be the unique one (up to an overall factor) by 
Reisenberger in j!04l . \^bc) defines the so-called Barrett-Crane intertwiner. Now the 
projector Pf nv in (JHUj) — the building block of the BF amplitude — can be written as 

P* nv = \Vbc) ($bc\ + orthogonal terms, (88) 

using the standard Gram-Schmidt construction of a basis in Inv [pi • • ■ pj. In other words 
Pf nv is the sum of 1-dimensional projector to the solutions of the constraints (J82|) plus the 
orthogonal complement. The solution to (J82)) is then unique (up to scaling) and can be 

28 In a detailed analysis Baez and Barrett ^2] justify this choice in a rigorous way. This restriction (im- 
posed by the so-called chirality constraint, Section l7T|l implies that the 'fake tetrahedron' configurations — 
corresponding to solutions of the constraints on the right of equation (JHSJ —are dropped from the state 
sum. 
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written as 




where the pij = jij ® j*j and we graphically represent bc) by 

The amplitude of an arbitrary simplicial complex is computed by putting together 
4-simplexes with consistent boundary connections and gluing them by integration over 
boundary data. If we do that we obtain the Barrett-Crane state sum on a fixed triangula- 
tion 

(bc) 

zbc(ja) = e n &f + x ) 2 n ^ n J\ L > w 

C : :\f\ -j. feJA eGjA \ / V\" \ / 



where we have made the replacement A p = (2j + l) 2 for p = j <g> j*. The vertex amplitude 
depends on the ten representations labeling the ten faces in a 4-simplex and it is referred 
to as lOj-symbol. Using the definition ()87j1 of the Barrett-Crane intertwiner, x ^ x , the 
lOj-symbol can be written explicitly in terms of 15j-symbols as 




where we represent the normalized S'0(4)-intertwiners in (|57jl in terms of the cor- 

responding pair of left-right S'£7(2)-intertwiners. 

A e denotes a possible edge amplitude which is not determined by our argument. The 
implementation of constraints in the path integral should be supplemented with the ap- 
propriate modification of the measure. This should affect the values of lower dimensional 
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simplexes such as face and edge amplitudes. Constraints (jH^j) act on each edge (tetra- 
hedron) separately; heuristically one would expect a Jacobian factor to modify and so 
determine the value of the edge amplitude A e . As we pointed out in Section lfTT| the reso- 
lution of this ambiguity is related to the question of how to define the correct anomaly-free 
measure. In Section El we will see how the GFT formulation provides a natural definition 
of A e . 

Constraints that involve different tetrahedra in a given 4-simplex — corresponding to 
(172)1 — are automatically satisfied as operator equations on the Barrett-Crane solutions. 
This can be checked using (j%2)) and (j7T)|) . 

7.4 The quantum tetrahedron in 4d and the Barrett-Crane in- 
tertwiner 

We have re-derived the Barrett-Crane model directly from a simplicial formulation of Ple- 
banski's action. The original definition makes use of the concept of the 'quantum tetrahe- 
dron' [THUS]. In that context, the analog of our £?'s is given by the bivectors associated 
with the triangles of a classical tetrahedron. A Hilbert space is defined using geometric 
quantization and the classical triangle bivectors are promoted to operators. This Hilbert 
space is reduced by implementation of constraints — quantization of the geometric con- 
straints satisfied by the classical tetrahedron — to the Hilbert space of the so-called 'quan- 
tum tetrahedron'. These constraints are precisely of the form 1)72)1 . For fixed triangle 
quantum numbers the state of the quantum tetrahedron is defined by the Barrett-Crane 
intertwiner ()87jl . There is however no systematic prescription for the construction of the 
state sum and even the single 4-simplex amplitude has to be given as an, although natural, 
ad hoc definition. 

As we have shown in the previous sections all this can be obtained from the path in- 
tegral approach applied to ()65j) . In addition to simplifying the derivation — no additional 
quantization principle is required — our framework provides a direct and systematic defini- 
tion of the state sum. In ^5) E] one can only access the constraint operator algebra of a 
single tetrahedron and conclude that this sub-algebra is closed. The full constraint algebra 
of (J82j) can only be studied in the formalism presented here. Finally, since our prescription 
is directly derived from an action principle it is conceivable that a rigorous derivation of 
the undetermined amplitudes (such as A e in (jnUjl ) could exist. These important questions 
remain to be investigated. 

7.5 An integral expression for the lOj-symbol 

In Section El we will present a discretization independent formulation of the Barrett-Crane 
model based on a GFT. The realization by Barrett that the vertex amplitude ()91)1 
admits an integral representation plays a key role in the construction. The integral formula 
is also important in the computation of the lOj-symbol asymptotics that we briefly describe 
in the next subsection. 
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The basic observation in the construction of an integral expression for ()91|) is that 
equation (|87|) has precisely the form (JHHj) (for n — 4) if we write the Spin(4) intertwiners 
as tensor products of SU(2) ones. Therefore, 



\Vbc) = J du 3l (u) <g> j 3 (u) <g> j 3 (u) ® j 4 («), (92) 

S*[/(2) 

where j'(m) denotes SU(2) representation matrices in the representation j. Each one of the 
five pairs of intertwiners in (|HTj) can be obtained as an integral ()92j) over SU(2). Each of the 
ten representation matrices j%k (i ^ k — 1- ■ - 5) appears in two integrals corresponding to 
the intertwiners at the node i and j respectively. Contracting the matrix indexes according 
to (}9"Tj) these two representation matrices combine into a trace Tr [jikiuiU^ 1 )] (ui G SU(2)). 
Parameterizing SU(2) with spherical coordinates on S 3 

Tr WW)] = ■= ( 2 ** + !)*i*(w, W), (93) 

where ^ is the azimuthal angle between the points yi,yk on the sphere corresponding to 
Ui and Uk respectively. We have also introduced the definition of the kernel Kj ik (yi,yk) in 
terms of which the Barret-Crane vertex amplitude (j91j) becomes 

5 

AvUik)= / \{dy l \{{2 3i k + l)K Hk {y l ,yk). (94) 



i<k 



Each of the five integration variables in S 3 can be regarded as a unit vector in R 4 . They 
are interpreted as unit normal vectors to the 3-dimensional hyperplanes spanned by the 
corresponding five tetrahedra. The angles ipik is defined by cos ipik = y%-yk and corresponds 
to the exterior angle between two hyperplanes (analogous to the dihedral angles of Regge 
calculus). These normals determine a 4-simplex in M 4 up to translations and scaling [105 . 



7.6 The asymptotics for the vertex amplitude 

The large spin behavior of the spin foam amplitudes provides information about the low 
'energy' or semi-classical limit of the model [T06J in the naive sense h — > oo while geometric 
quantities such as the area (J7J) are held fixed. Evidence showing a connection between the 
asymptotics of the Barrett-Crane vertex and the action of general relativity was found by 
Crane and Yetter in jl()7j . 

A computation of the asymptotic (large j) expression of the Barrett-Crane vertex am- 
plitude for non-degenerate configurations was obtained by Barrett and Williams in [108 . 
They computed A v (jn-) for large jik by looking at the stationary phase approximation of 
the oscillatory integral (|94j). The large spin behavior of the vertex amplitude is given by 



7T 



S Regge ((T) + K- 



D 



(95) 
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where the sum is over geometric 4-simplexes cr whose face areas are fixed by the spins. The 
action in the argument of the cosine corresponds to Regge action which in four dimensions 
is defined by S Regge (a) = J2i<kAik V'ifc( cr ) where is the area of the i/c-triangle. P(a) 
is a normalization factor which does not oscillate with the spins. D is the contribution 
of degenerate configurations, i.e. those for which some of the hyperplane normals defined 
above coincide. However, in a recent paper |lU9j Baez, Christensen and Egan show that 
the term D is in fact dominant in the previous equation, i.e. the leading order terms are 
contained in the set of degenerate configurations! This has been later confirmed by the 
results of Freidel and Louapre and Barrett and Steele 

7.7 Area and Volume in the Barrett-Crane model 

Using the representation of the discrete B of Section 1731 we can define the operator corre- 
sponding to the square of B t corresponding to any triangle t 6 A, namely 

B tIJ B{ 3 = -X u (U f )X IJ (U f ), 

where / is the face dual to t G A. Acting on the modes of the expansion this is simply 
the Casimir 

B tIJ Bl J = [/(/ + 1) + f{f + 1)] 1. (96) 

It is easy to check that this operator commutes with all the constraints and therefore is 
well defined on the space of solutions of the constraints. Its geometrical meaning is clear 
if we recall that at the classical level B t represent the bivector associate to the triangle 
t e A. Once the constraints are imposed (|96|) is proportional to square of the area of the 
corresponding triangle. The simplicity constraint implies je = j* = j so that the area of 
the triangles have discrete eigenvalues a,j given by 

aj cc x/iO' + l), (97) 

in agreement with the result of quantum geometry (J7J) 29 . 

Can we define the analog of volume operator of quantum geometry? The candidate for 
the square of such operator are 

where Ji, J2, and J 3 are the operators defined in corresponding to three different 
triangles in the tetrahedron. The operator U + vanishes identically on the solutions of the 
quantum constraints. In fact it can be expressed as the commutator of the constraints 
([82)1 . This is the chirality constraint of [To] . 

One would like to define the volume operator as the square root of £/_; however, Z7_ 
is not a well defined operator on the Hilbert space of the quantum tetrahedron since it 

29 Given a triangulation A of M. and the induced triangulation As of a slice S C M. any spin foam 
denned on J\ induces a spin network state on a graph j s dual to As- Links of 7^ are dual to triangles 
in As . These triangles play the role of the surface S in the definition of the area operator @ . 
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does not commute with the simplicity constraints. In other words the action of the volume 
operators map states out of the solution of the constraints. Only when the dimension 
of Inv [ji (gi j'2 <8> J3 <8> ji] is one the commutator vanishes (this is in fact a necessary and 
sufficient condition). The volume operator is well defined in this subspace but it vanishes. 
In conclusion, generically the volume operator is not defined in the Barrett-Crane model. 

7.8 Lorentzian generalization 

A Lorentzian generalization of the Riemannian Barrett-Crane vertex amplitude was pro- 
posed in jTQJ - Using the spin foam model GFT duality of Section |H1 a generalization of 
the full model (including face and edge amplitudes) was found in [13J. The relevant repre- 
sentations of the Lorentz group are the unitary ones. Unitary irreducible representations 
of SL(2, C) are infinite dimensional and labeled by n G N and p G M + . The simplicity 
constraints select the representations for which n = 0. The triangle area spectrum is given 



There is still a minimum eigenvalue but the spectrum is continuous. The derivation of the 
model can be obtained from £0(3, 1) Plebanski's action by generalizing the computations 
of Section We will study this model in Section 19.31 were we show that the Lorentzian 
extension is not unique. In fact a new Lorentzian model can be defined following the 
Barrett-Crane prescription |112j . 

7.9 Positivity of spin foam amplitudes 

Baez and Christensen |113j showed that spin foam amplitudes of the (Riemannian) Barrett- 
Crane model are positive for any closed 2-complex. For open spin foams they are real and 
(if not zero) their sign is given by (— 1) 2J where J is the sum of the spin labels of the edges of 
the boundary spin networks. This is a rather puzzling property since (as the authors point 
out) this seems to imply the absence of quantum interference. Positivity of the Lorentzian 
model seems to hold according to numerical evaluations of the vertex amplitude. The 
positivity results have been generalized to generic 2-complexes (not necessarily dual to 
a triangulation) by Pfeiffer in j!14j . who also showed that the sign (— 1) 2J is not really 
present if one chooses the correct bases of intertwiners. 

7.10 Degenerate Sectors of Spin(4) Plebanski's formulation 

In jSHj Reisenberger showed that, solving the constraints (|66|) in the degenerated sector — 
i.e., where e = ^opqrB?^ 'B^e^ up(T = — and substituting the solution back into the 
Plebanski action, the theory reduces to two sectors described by the actions 



by 



a p oc 



(99) 




(100) 



48 



where the upper index r (respectively £) denotes the self-dual (respectively anti-self-dual) 
part of B and A in the internal space, and V G SO (3). Let us concentrate on the sector 
with the minus sign in the previous expression. Then it is straightforward to define the 
discretized path integral along the same lines as BF theory in Section [7"T1 The result is 

Z(A) = / II dBf ] dv f J] dgidgl /MW/^V] , (101) 

where dg l e1 dg r e1 and dvf are defined in terms of the SU{2) Haar measure. Integrating over 
the B field we obtain 

z(A) = f n rf^i n cfo, • • ■ s e >/G4 • • ■ gixwi ( io2 ) 

^ egA* /eA* 

where the delta function 5^ denotes an SU(2) distribution. 

Integration over g l e) g r e and Vf yields a spin foam model where only face representations 
are constrained to be simple while intertwiners are arbitrary. Explicitly 



1: 



Zdeg (A)= £ En" 

C f :{f}^j f C e :{e}^{u e } veA* 




(103) 



where pj = jf <E> jj and jf G Irrep[S'f/(2)]. This is precisely the spin foam obtained in [TTj ! 
This model was obtained as a natural modification of the GFT that defines a variant of 
the BC model. Here we have rediscovered the model from the simplicial quantization of 
Sfe defined in ([100)1 . This establishes the relation of the model with a classical action! It 
corresponds to spin foam quantization of the '— ' degenerate sector of 5*0(4) Plebanski's 
theory. 

The + sector action ()100|) can be treated in a similar way. The only modification is 
that of the subgroup. Instead of using the diagonal insertion defined above one has to 
define u G SU(2) C Spin(4) so that ug = (ug i ,u~ 1 g r ). This selects representations of the 
form p = j (gi j instead of p = j <8> j* for faces. 

Notice that the allowed 4-simplex configurations of the model of Section 17.101 are fully 
contained in the set of 4-simplex configurations of the model obtained here. 

8 The group field theory (GFT) ansatz 

In this section we present the motivation and main results of ^Sl El where the duality 
relation between spin foam models and group field theories is established. This is the 
formulation we referred to in Section 16.21 as one of the possible discretization independent 
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definitions of spin foam models. The result is based on the generalization of matrix models 
introduced by Boulatov 86] and Ooguri jHTJ dual to BF theory in three and four dimensions 
respectively. This formulation was first proposed for the Riemannian Barrett-Crane model 
in ^Tj and then generalized to a wide class of spin foam models |16| I17j. 

Given a spin foam model defined on an arbitrary 2-complex J7a (dual to a triangulation 
A) -thus the partition function Z\J^\ is of the form (}2T]) - there exists a GFT such that the 
perturbative expansion of the field theory partition function generalizes ()21|) to a sum over 
2-complexes represented by the Feynman diagrams J of the field theory. These diagrams 
look 'locally' as dual to triangulations (vertices are 5-valent, edges are 4-valent) but they 
are no longer tied to any manifold structure [881 • 

We now motivate the duality using what we know of the model defined on a fixed 
simplicial decomposition. The action of the GFT is of the form 

J[0] = / o [0] + ^[0], (104) 

where Io[4>] is the 'kinetic' term quadratic in the field and V[<p] denotes the interaction 
term. The field <fi is defined below. The expansion in A of the partition function takes the 
form 

r \ N 
z = J me -m = E z[Jn] , {m) 

where Jn is a Feynman diagram (2-complex) with N vertices and Z\J^\ is the one given 
by (|2ip. In fact, the interaction term V[(f>] is fixed uniquely by the 4-simplex amplitude of 
the simplicial model while the kinetic term Io[4>] is trivial as we argue below. 

Expression ()105j) is taken as the discretization (and manifold) independent definition 
of the model. Transition amplitudes between spin network states on boundary graphs 71 
and 72 are shown to be given by the correlation functions 

/ e- m , (106) 
71 72 

where boundary graphs are determined by the arrangement of fields in the product. The 
field can be defined as operators creating four valent nodes of spin networks |115t l§2~] . In 
[H2] correlation functions (0 • • • 0) are interpreted as a complete family of gauge invariant 
observables for quantum gravity. They encode (in principle) the physical content of the 
theory and can be used to reconstruct the physical Hilbert space in a way that mimics 
Wightman's procedure for standard QFT \IW\ I117j . 

How do we construct the GFT out of the spin foam model on a fixed 2-complex? In 
f!105|) the combinatoric of diagrams is completely fixed by the form of the action ()104|) . 
To construct the action of the GFT one starts from the spin foam model defined on a 
2-complex J7a dual to a simplicial complex A. The sum over 2-complexes in ()106|) contains 
only those 2-complexes that locally look like the dual of a simplicial decomposition. 

Spin foams on such 2-complexes have edges e with which are associated a tensor product 
of 4 representations p\ <g> pi <g> ps ® p±, since they bound four colored faces in J7a- If we 



50 



want to think of the edge e as associated to the propagator of a field theory then such a 
propagator should be a map 



V '■ Pl <£> 92 <8> P3 <8> P4 — ► Pi <£> 92 <8> P3 ® P4- 



(107) 



According to Peter- Weyl theorem, elements of pi (8> ■ • • <8> P4 naturally appear in the mode 
expansion of a function (f)(xx, ■ ■ ■ ,x^) for X{ G G. Spin labels arise as 'momentum' vari- 
ables in the field theory. Intertwiners assigned to edges in spin foams impose compatibility 
conditions on the representations. In the context of the GFT this is interpreted as 'momen- 
tum conservation' which is guaranteed by the requirement that the field be 'translational 
invariant'. 

(f)(xi,x 2 ,x 3 ,x i ) = (f)(x 1 g,x 2 g,x 3 g,x 4 g) V geG. (108) 

One also requires <fi to be invariant under permutations of its arguments. One can equiva- 
lently take the field (f> to be arbitrary and impose translation invariance by acting with P 
defined as 

P(j)(x 1 ,X2,x 3 ,x 4 ) = / dg (j)(x 1 g,x 2 g,x 3 g,x A g). (109) 



All the information of local spin foams is in the vertex amplitude or fundamental atom 
(4-simplex) amplitude of Section ET2l Therefore, the kinetic term Io[<p] is simply given by 



^ J d 4 x [P0(xi,x 2 ,x 3 ,x 4 ) 



(110) 



Since vertices are 5-valent in our discretization the interaction term should contain the 
product of five field operators. This is a function of 20 group elements. If we use the 
compact notation := (j){x\, • • • , £4), then the general 'translation invariant' form is 



v[4>] 



d x V(xij[x 



P0(Xii)P0(x2j)P0(x 3i )P0(x4i)P0(a;5i 



(111) 



where V is a function of 10 variables evaluated on the 'translation invariant' combinations 
atij := Xijlxji} -1 . For local spin foams the function V{otij) is in one-to-one correspon- 
dence with the fundamental 4-simplex (atom) amplitude (|19|) . If we represent the field 

0(xi, X21 xs, £4) by <^f^ we can write the action ()104|) . using (jllOjl and (|111|) . as 

P4> 




S[<f>] 













1 d 4 x 1 






3 















d 20 x 
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where the straight lines represent the field arguments (in the case of the interaction term, 
the 20 corresponding x^- for i ^ j = 1, • • • , 5). 

The precise way in which 2-complexes are generated as Feynman diagrams of the GFT 
will be illustrated in the following section. We have shown that any (local in the sense 
of Section lr>.2j) spin foam model can be given a discretization independent formulation in 
terms of a uniquely determined GFT theory. 

9 The Barrett-Crane model and its dual GFT-formulation 

In addition to providing a discretization independent formulation of the Barrett-Crane 
model, the GFT formulation provides a natural completion of the definition of the model. 
In Section [7| (equation (J§U|) ) we noticed that the model on a fixed discretization is defined 
up to lower dimensional simplex amplitudes such as that for faces / (dual to triangles) and 
edges e (dual to tetrahedra). The GFT formulation presented here resolves this ambiguity 
in a natural way. This normalization of the Barrett-Crane model was also obtained in |118j 
using similar techniques but on a fixed triangulation. 



In this section we introduce the general GFT action that can be specialized to define the 
various spin foam models described in the rest of the paper. 

9.1.1 The field theory action and its regularization 

Consider the Lie group G corresponding to either Spin(4) or SL(2, C) - for the GFT dual 
to Riemannian or Lorentzian Barrett-Crane model respectively. The field <f>(xi, x%, X3, 24) 
is denoted 4>(xi) where i = 1 ... 4, and x,- t e G, symmetric under permutations of its 
arguments, i.e., <p(xx, x 2 , x 3 , x±) = <p(x a m,x a f 2 ),x a f^,x a m), for a any permutation of four 
elements. Define the projectors P and R as 



where U C G is a fixed subgroup, and dg and du are the corresponding invariant measures. 
The projector P imposes the translation invariance property (|108|) . 

Different choices of the subgroup yield different interesting GFT's. When U — {1} 
(R = 1) the GFT is dual to BF theory and we get the model of [HZJ - The GFT is dual 
to the Riemannian BC model for U = 577(2) C Spin(A). Similarly for the GFT dual to 
the Lorentzian models the subgroups U C SL(2, C) are U = SU(2) (leaving invariant a 
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(113) 



G 



and 




(114) 
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time-like direction) or U — SU(1, 1) x Z2 (for an invariant space-like direction); they result 
in two different models. 

The GFT action is is of the general form ()104j) and simply given by 



S[</>] = J dx, t [P0(x,)] 2 + A J JJPRP^(xy), (115) 

where i,j = 1...5,i^j and dx is an invariant measure in G (the normalized Haar 
measure in G is compact). 

Strictly speaking the operators P and R are projectors only when the corresponding 
groups are compact. Formally we have 

P 2 = vol G x P, and R 2 = vol v x R, (116) 

where voIq, and volu denote the volume of G and U, respectively. These volume factors 
can be taken to be one when the G and/or U are compact by using Haar measures. When 
G and/or U are non compact the factors are infinite. This is a rather simple technical 
problem with which we shall deal later. Essentially one must drop redundant projectors 
in the functional integral. 

The partition function can be computed as a perturbative expansion in Feynman dia- 
grams Jjsii namely 

r \ N 

Z = / = TT^MJn), (H7) 

J sym(J N ) 

where N is the number of vertexes in Jn and sym(j7}v) is t ne standard symmetry factor. 

In what follows we explain the structure of the Feynman diagrams of the theory. If we 
use the notation §{x\j) = <fi(xi2, #13, £14, £15) we can write the action as in (jllUj) and (jlllj) 

S[(f>] = - J dxidx\ (f)(xi) K,(xi, x'j) 4>{x'j) 
A 



+ 5! / dXij 0( Xl ^ ^ x<2 ^ ^( x ^j) ^( x 4j) <j>(x5j), (118) 

where i 7^ j and the kinetic /C(xi,x-) and interaction V{xy) operators a explicitly given 
below. The kinetic operator /C(xj,x-) is given by 

K,{x h x\) = J d 9 d 9" Yl 5 ( Xig " 9 '~ lx 'v(i)} 

a G2 8=1 

= zJ / d9 ' d9 Tl S ( X i9 X 'a(i)) 

a G2 i=l 

= vol G xY^ J dg H6(x i9 x>-i), (119) 



G 
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where the g" and g' integrations correspond to the action of the projectors P in ()115j) . 
Redefining the integration variables g = g" g'~ x we obtain the second line in the previous 
equation. The vo\q factor comes from the g' integration as in (|116|) . We regularize the 
kinetic operator by simply dropping one of the G-integrations in the previous expression, 
namely 

r 4 

JC(x t ,x[) =J2 J dg n^<»)- ( 12 °) 

a q i=l 

P acts by projecting the field into its 'translation invariant' part P0(xj). The action 
()115j) depends only on the gauge invariant part of the field, namely S[<f>] = S[P<p] (recall 
1)108)) ). The inverse of K (in the subspace of right invariant fields) corresponds to itself; 
therefore the propagator of the theory is simply 

V(x i ,x' i )=JC(x i ,x , i ). (121) 

The propagator is defined by 4 delta functions (plus the symmetrization and the integration 
over the group) and it can be represented as shown on the right diagram of Fig. [7| The 
potential term ()118|) can be written as 

V{xij) = dgidg-duij Y[5(xjlg' i u ij g^ 1 g j u j ig'r 1 x ij ), (122) 

' J i<j 

where the Uij G U correspond to the action of R in ()115|) and g^g^ G G to the action of 
the corresponding P's, respectively. It is easy to check that in the evaluation of a closed 
Feynman diagram the g" s can be absorbed by redefinition of the g^s in the corresponding 
adjacent propagators. In this process, the g"s drop out of the integrand and each integral 
over g\ gives a voIq factor. This means that the second P projector in (|115|) is redundant 
in the computation of ()117|) . The regularization is analogous to the one implemented in 
equation ([120)1 : we drop redundant P's. The regularized vertex amplitude (for a vertex in 
the bulk of a diagram) is then defined as 

V(xij-) = — j dgiduij Y^Sixjluyg^gjUjiXij). (123) 
' J i<j 

In the case of open diagrams g"s remain at external legs. 

There is still a redundant gi integration in (j 123)1 which introduces a potentially infinite 
v o\q factor in the vertex amplitude. Notice that the ([123)1 depends on the 'translational 
invariant' combinations g^gj so that one of the (7-integrations is redundant. In other 
words, one can absorb one of the by redefining the remaining four using the invariance 
of the measure. In the non compact case this would yield another vole infinite factor. The 
regularization now consists of removing an arbitrary <?j from the expression of the vertex 
amplitude. This results in the regularization scheme proposed by Barrett and Crane in 
[TU] . The regularization presented here can be applied to any non compact group model on 
a lattice and can be regarded as a gauge fixing condition for the internal gauge invariance 
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Figure 7: The structure of the propagator and the interaction vertex. Each line represents 
a delta function as the ones in the integrand of (|120|) in the case of the propagator and 
(|123J) in the case of the vertex. The shaded box in the propagator represents the sum 
over permutations a in (|120|) . This sum is over diagrams including crossings, we have 
represented three of this terms in the diagram below the propagator. 

(recall Section 16.1)1 . For notational simplicity we do not implement the regularization 
explicitly in ()123j) . The structure of the vertex is represented on the left diagram in Fig.0 
Each of the lines represent a delta function appearing in the integrand of (j!23)) (compare 
with[TT2J). 

The Feynman diagrams of the theory are obtained by connecting the 5-valent vertexes 
with propagators (see Fig. EJ). At the open ends of propagators and vertexes there are the 
four group variables x corresponding to the arguments of the field (in addition to the inte- 
gration variables gi and u%j coming from P and R respectively). For a fixed permutation a 
in each propagator, one can follow the sequence of delta functions with common arguments 
across vertexes and propagators. On a closed graph, each such sequence must close. By 
associating a surface to each such sequence of propagators, we construct a 2-complex J 
[TT] . Thus, by expanding in Feynman diagrams and in the sum over permutations in (|119)) . 
we obtain a sum over 2-complexes. Each 2-complex is given by a certain vertex-propagator 
topology plus a fixed choice of a permutation on each propagator. 

9.1.2 Evaluation in configuration space 

Combining (jl21j) and (j!23|) and integrating over internal variables X{j the sequence of delta 
functions associated to a face / G J reduces to a single delta function. Denoting v\ . . . v n 
the ordered set of n vertices bounding /, and the edge connecting Vi with Vj, the delta 
function corresponding to the face / becomes 

where g eij G G is the integration variable in (jl2U|) associated to the edge e^, g Vi . G G 
denotes the integration variable in (j!23)) associated with the leg ij of the vertex Vi connected 
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to the edge e^, and u f Vi . E U is the corresponding subgroup integration in one of the delta 
functions in (|123|) . We use the supra index / to emphasize that the it's appear on single 
delta functions in (|123J) and therefore do not contribute to more than one face. In the 
previous equation, the product of group elements between brackets correspond to the 
vertex contribution to the face (see ()123j) ) while the g ei . comes from the corresponding 
propagators ()121j) . also we have that g eij = g~\ Using (j!24|) the amplitude A(J) of an 
arbitrary closed 2-complex becomes a multiple integral of the form 

A{J) = j \{dg e \{dg v dg v (125) 

e v 

]Jduidui6([ui ln g-lg n y v Jg e Jui 2 ^^^ 
f 

Notice that if the subgroup U = {1} the previous expression coincides with the BF am- 
plitude (jfil)) if we make the change of variables g' e .. = gv^gei-g^- For U = {1} we recover 
Ooguri's GFT [HZ]. 



9.1.3 The spin foam representation 

The spin foam representation for the amplitude can be obtained by expanding the delta 
functions in terms of irreducible unitary representations of G, namely 

5(g) = J2\Tr[p(g)}, (126) 
p 

where p labels unitary irreducible representations and the rest of the notation is that of 
()31|) when G is compact. In the non compact case representations are infinite-dimensional 
so a formally equivalent expression holds where A p correspond to the so-called Pancharel 
measure |112j . 

Using ()126|) in ()124j) the face contribution is 



J2 A p Tr M<J • P(9vi n 9vJ ■ P« 2 ) ■ P{9eJ ■ P«J ■ P{9v2l9vJ- 
p 

P« 3 ) • P(9 e j ■ ■ •/>«(„_!)) • P(9vX-i)9vJ ■ P( U U ■ P(9ej] , (127) 

where we have dropped the integration symbols for simplicity. The it's appear just once 
per face so we can perform the w-integrations independently of other faces. Notice that 
j dup(u) is the projection onto the invariant subspace under the action of U in the Hilbert 
space TC P . When U = {1} the projector is the identity. In the other cases the subspace 
turns out to be 1-dimensional. Consequently the projector can be written as 

du p{u) = \w p ) (w p \ , (128) 



u 
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where \w p ) is the corresponding invariant vector. Equation (| 12 7(1 becomes 

Yl A p ( w P \p(9vL9v 12 ) w p) ( w p\p(9 ei2 ) w p) ( w p\ P(9v 2 ]9v 23 ) w p )--- (w p \ p(g enl ) w p ) . (129) 
p 

There are two types of terms in the previous equation: those corresponding to edges 
for which the representation p is evaluated at the single g eij group variable, and those 
corresponding to vertices for which p is evaluated at the product g~^g Vik . Let us integrate 
over the edge variables first. The group element g e associated to an edge e G J appears 
four times as there are four delta functions in the propagator (j!21j) . Therefore, integrals 
over such g's have the general form 



MpI, ■ ■ ■ ,pi) = I (w p i\p 1 (g) w p i) (w p 2\p 2 (g) w p 2) (w p 3\p 3 {g) w p s) (w p 4\p 4 (g) w pi ) , 

(130) 



G 



and define the edge amplitude. If we define the kernel K p (g) as 

K p (g) = (w p \p{g)w p ), (131) 
then using the subgroup invariance of the |w p )'s the previous equation becomes 

A e ( Pl , . . . ,p 4 ) = volu x J dy K pi (y)K p2 (y)K ps (y)K pi (y), (132) 

G/U 

where the integration is over the homogeneous space G/U. Integration over the group 
elements associated to vertices are of the general form 



5 

A v (p ik ) = / \\dyt \\K Pik {y u y k ). 



4 i=l i<k 



(G/U) 

This corresponds to the vertex amplitude. When U is non compact the preceding expres- 
sions have to be regularized as in Subsection 19 . 1 . II bv dropping an arbitrary ^-integration. 
The amplitude of a 2-complex J is then given by 

A ^= E U\ f llA(pi,---,pdl[A v (p 1 ,...,p w ) (133) 

Cf-{f}-*{p} fe J e£j v£j 

As we shall see in the following the presence of R in (|115|) can be interpreted as imposing 
the Barrett-Crane constraints on the GFT dual to BF theory. 

9.1.4 Boundaries 

The amplitude of an open diagram, that is, a diagram with a boundary, is a function of 
the variables on the boundary, as for conventional QFT Feynman diagrams. The boundary 
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of a 2-complex is given by a graph. To start with, the amplitude of the open diagram is 
a function of 4 group arguments per each external leg. However, consider a surface of an 
open 2-complex and the link ah of the boundary graph that bounds it. Let a and b be the 
nodes on the boundary graph that bounds ab. The surface determines a sequence of delta 
functions that starts with one of the group elements in a, say g a , and ends with one of the 
group elements in b, say g b . By integrating internal variables, all these delta functions can 
be contracted to a single one of the form S(g a ■ ■ ■ g^ 1 ) = 6(g^g a • ■ ■ ). We can thus define 
the group variable p ab = g b x g a -, naturally associated to the link ab, and conclude that the 
amplitude of an open 2-complex is a function A(p ab ) of one group element per each link of 
its boundary graph. A(p ab ) is gauge invariant as can be easily checked, and p a b represents 
the discrete boundary connection in the language of Section 13.21 

In "momentum space" , boundary degrees of freedom are encoded in spin-network states. 
That is, if s is a spin network given by a coloring of the boundary graph, then 

A(s) = J dp ab i) s { Pab ) A(p ab ), (134) 
where ip s (Pab) is the spin network function [2111211 • The formula can be inverted 

A( Pab ) = J2 A (s)Mpab)i (135) 

s 

where the sum is performed over all spin network states defined on the given boundary 
graph. The argument presented here applies in general for any GFT of the kind introduced 
in Section |HJ In the case of the GFT's defined here p ab = g b x g a will appear between U- 
integrations. So strictly speaking, for each p ab = g b l g a there correspond two point y a and 
y b in the homogeneous space G/U and one cannot say that the amplitude depend on the 
G-connection. 



9.2 GFT dual to the Riemannian Barret-Crane model 

The representation theory of Spin{4) is particularly simple due to the fact that Spin(A) = 
SU (2) x 577(2). Unitary irreducible representations of p of Spin(4) are labeled by two half- 
integers ji and j r , and are given by the tensor product of unitary irreducible representations 
of SU(2), namely 

Phir=3i®3r- (136) 
In terms of representation matrices we have 

R jljr (g)mm' qq > = DtMD^igi), (137) 

where g = (gi,g r ) £ Spin{£), and -D^ m , are S'?7(2)-representation matrices. The subgroup 
U of previous section is taken as U — SU(2) C Spin(A) defined by the diagonal action on 
Spin(A) as follows 

gu = (giu,g r u) for u E U. (138) 
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The projector ()128|) can be explicitly evaluated 



v 



du p jljr (u) = J du D^ m ,(u)D j qq ,(u) = ^-j5 mq 5 m/q , = \w n ) (w n | , (139) 



where we have used the ort honor mality of 577(2) representation matrices (Footnote irU|). 
The previous equation confirms that the invariant subspace for each p^ T is 1-dimensional 
as anticipated above. An immediate consequence of (|139|) is that the kernel K p in (J131j) 
becomes 

Kj{g) = {w 3 \p jl0 {g) Wj ) = rJ^Tr [D^g l9r )] . (140) 

The factor 8j>j restricts the representation to the simple representations ji = j* = j. 

The projector R in f)115j) imposes the Barrett-Crane simplicity constraints of equation 
(|83jl ! . Simple representations appear in the harmonic analysis of functions on the 3-sphere, 
53. The fact that we encounter them here is not surprising since 53 is the homogeneous 
space 5 3 = Spin(4)/SU(2) under the SU{2) diagonal insertion (|138|) . The presence of R 
in f)l 15|) projects out those modes that are not simple or spherical. Indeed, (w^ Pj>j*(g) Wj) 
can be thought of as a function on 53: (w } \ pj>j* (g) Wj) depends on the product gig r G SU (2) 
which is isomorphic to 53. We will see below that in fact R imposes all the constraints of 
Section fixing in addition the value of the so far undetermined edge amplitude A e . 

The invariant measure on SU(2) can be written as a measure on 53 induced by the 
isomorphism SU(2) — > S3. We can parameterize h £ SU(2) as h = y^ h) r^ where = iak 
for k — 1, 2, 3, {(Tfe} is the set of Pauli matrices, and r = 1. h E SU (2) implies y^y^n = 1 
(indexes are lowered and raised with 5%), i.e., y £ S3. In this parameterization the SU(2) 
Haar measure becomes 

dh^dy = ^dif6(y*y lt -l), (141) 

We can simplify the measure using spherical coordinates for which an arbitrary point y G 5 3 
is written 

y = (cos(^), sin(?/')sin(6 l )cos((/)), sin('0)sin(6 l )sin(0), sin(^)cos(6 1 )) , (142) 
where 0<^<7r, 0<#<7r, and < <fi < 2ix. The Haar measure becomes 

dy = At sin 2 (?/>) sin(0) dip dd dtp. (143) 
2n 2 

Using the known formula for the trace of SU(2) representation matrices it is easy to check 
that ()140j) becomes 

Wy) = Si ' j (144) 

3WvJ 2j + 1 sin(^) ' 1 ; 

where iph is the value of the coordinate ip in (|142|) corresponding to y G SU(2). Notice 
that we have rediscovered the kernel (J93|) which together with equation (J94|) establishes 
the claimed duality between the Barrett-Crane model and the GFT of this section. 
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Now, the general edge and vertex amplitudes f ()132j) and (|133p ) reduce to multiple 
integrals over S3. These can be interpreted as the evaluation of diagrams on the sphere 
with propagator given by f|144jl . This evaluation is referred to as relativistic spin network 
evaluation in the literature. In the case of the edge amplitude the evaluation corresponds 
to a 84 graph with four links labeled by the four spins ji, ■ ■ ■ , J4. In the case of the vertex 
amplitude the evaluation corresponds to the 75 spin network of Figure 0] for which the ten 
links are labeled by the corresponding ten spins ji, . . . ,jiQ. There is no intertwiner label 
at nodes as the projection R fixes the intertwiners to the unique Barrett-Crane one. 

The amplitude A e (ji, . . . , j'4) can be computed using (jl4()|) and (jHHJ), it follows that 

Ae{ji,...,j4) = -r. r. . . . 1 , 145 

dim [j 1 ® j 2 ® 33 ® 34 

where dim [Inv(ji ® j 2 ® 33 <S> ji)] denotes the dimension of the trivial component of Tt^ <E> 
TCj 2 <8> Tlj 3 ® TCj 4 and dim [ji ® j 2 <8> J3 <8> J4] the dimension of the full space 30 . So finally, 

= e n^/ + x ) 2 n ^(^' • • • >ti n ^o - i. ■ ■ ^ 

c f -.{f}->{j} fej eej vej 

where the vertex amplitude is given by (|133J) using (jl44j) . This is the Barrett-Crane model 
([90)1 where the values of A e and Af have been fixed by the GFT formulation. 

9.2.1 Finiteness 

In this section we prove that the infinite spin sum in (|146j) is convergent if the 2-complex J 
is finite and non-degenerate. In order to do this we need to construct bounds for the terms 
in the sum (edge and vertex amplitudes). These bounds are provided by the following 
lemmas. 

Lemma 9.1. [18] For any subset of k elements j\...j K out of the corresponding four 
representations appearing in A e (ji, . . . ,j±), the following bounds hold: 

14/. . \ 1 . 1 , fl for k < 3 

\A e {j x ,...,j A )\ < , where a K = I 3 . 

n(2j, + l)- 1 1 f0 

i=l 

The inequality for k < 3 is sharp. 

The proof of this lemma is elementary and can be found in ^H] ■ 

Lemma 9.2. 1181/ For any 4 spins j\ . . .j± labeling links converging at the same node in 
the relativistic spin-network corresponding to A v (j\, . . . , jio) the following bounds hold: 

1 



\A v (ji, . . -Jio)\ < 



(2j 1 + l)(2j 2 + l)(2j 3 + l)(2j4 + l) 



30 There was a mistake in the computation of the edge amplitude in due to the propagation of a 
typo in The erroneous expression of the edge amplitude contained dim \jx ® j'2 ® ja ® ji] in the 

denominator. 
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from which follows 



\A v {ji, ■ ■ ■ ,jw)\ < 



((2 Jl + l)...(2j 10 + l)) 2 / 5 ' 



Definition 9.1. A 2- complex J is said to be degenerate if it contains some faces bounded 
by only one or two edges. 

Theorem 9.1. L 18l The state sum A(J~), ftl4b\ ), converges for any non degenerate 2- 
complex J . 



Proof. The amplitude ()146j) can be bounded in the following way: 

\mj)\ < n e + i)) 2_|n/ ^ n/ = n e + i )) 2 " fn/ > ( i4? ) 

feJ j f feJ if 

where rif denotes the number of edges of the corresponding face, and we have used the fact 
the number of edges equals the number of vertices in a face of J . The term (2j + 1) 2 in (11471) 
comes from the face amplitude, (2j + l) - *™^ from Lemma f9. II (k, = 4), and (2j + l)~io n / 
from Lemma ()9.2j) k = 10. Notice that if the 2-complex contains faces with more than 
two edges the previous bound for the amplitude is finite, since the sum on the RHS of the 
previous equation converges for nj > 2. □ 

The fact that the sum over spins converges for a fixed (non degenerate) 2-complex is an 
very encouraging result. It means that for a fixed 2-complex the sum over the infinite set 
of 4-geometry configurations, represented by the corresponding spin foams, is well defined. 
Even though the result was obtained in the context of the GFT formulation, it can be 
applied to the model defined on a simplicial decomposition of the space-time manifold. This 
is because f!146p refers to the amplitude of a single Feynman diagram. In this context, the 
result is very encouraging and appears as a first necessary step for the construction of the 
discretization independent generalized projector and physical scalar product by studying 
the refinement limit of Section lfi~!fl On the other hand, if the discretization independent 
formulation is to be obtained directly by means of the GFT formulation, this result means 
that the perturbative expansion in A ()105|) is finite order by order. The previous statement 
has not been completely proved as degenerate 2-complexes — for which our proof does not 
hold — will appear in the A-expansion. However, we expect these diagrams to be finite also. 
A proof of this is lacking at the moment but it should be feasible by strengthening the 
bounds we have used in this section. 



9.3 Lorentzian models 

Unitary irreducible representations of SL(2, C) appearing in the general expression ()126|) 
are the ones in the so-called principal series. They are labeled by a natural number n 
and a positive real number p. Unitary irreducible representations of SL(2, C) are infinite 
dimensional. Those in the principal series are defined by their action on the linear space 
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V np of homogeneous functions of degree (^y 2 — 1, ~ n j lp — 1) of two complex variables z\ 
and z%. Due to the homogeneity properties of the elements V np they can be characterized 
by giving their value on the sphere \z\\ 2 + \z2\ 2 = 1. The latter is isomorphic to SU(2). 
The so called canonical basis is defined in terms of functions on SU (2) and is well suited 
for studying the following model. The relevant facts about SL(2, C) representation theory 
and the notation used in this section can be found in the appendix of reference |112j (for 
all about SL(2, C) representation theory see |119j ). 

9.3.1 GFT dual to the Lorentzian Barrett-Crane model 
Equation ()128j) becomes, in this case, 



du p np {u) = J duD^ e , m ,(u) = 5 n0 5 m 5 m = \w np ) (w np \ , (148) 

U 5(7(2) 

in terms of the canonical basis |112j . The SU(2) invariant vectors \w np ) are given by 
(21,22 \w np ) = 5 n o(\zi\ 2 + \z2\ 2 )% p ~ 1 which is indeed a homogeneous function of degree 
(f-l,f-l). 

An immediate consequence of ()148|) is that the kernel K p in (j!31|) becomes 

Knp(g) = (w np \ p np (g) w np ) = 5 n0 D^ 00 (u). (149) 

The kernel (w 0p \ po p (g) w 0p ) can be thought of as a function on H + = SL(2, C)/SU(2) 
(the upper sheet of 2-sheeted hyperboloid in Minkowski space). An arbitrary point y G H + 
can be written in hyperbolic coordinates as 

y = (cosh(?7), sinh(^)sin(6 l )cos((/)), sinh(^)sin(6 l )sin(0), sinh(^)cos(6')) , (150) 

where 0<r7<oo,0<#<7r, and < <p < 2tt. The invariant measure in these coordinates 
is 

dy = —— sinh 2 (r/) sin(^) dr]d6d(f). (151) 
27H 

In these coordinates, (11491) becomes 



K M = (152) 
p smh(?7j / ) 

where t] y is the value of the coordinate r\ corresponding to y G H + (this the generalization 
of (J9"3j) proposed in [ID])- Finally, the amplitude of an arbitrary diagram ()133J) becomes 

A{J)= J H P }dp f l[A e ({p e })l[A v ({p v }), (153) 



where the formal sum in the general expression (j!33)l becomes a multiple integral over 
the coloring Cf : {/} — > {p} of faces, and the weight A p now is given by the Pancherel 
measure, p 2 dp, of SL(2, C). The vertex amplitude was proposed in [TO], the full state sum 
amplitude including the previous normalization was obtained in 
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9.3.2 Finiteness 



The amplitude of a non degenerate 2-complex turns out to be finite as in the Euclidean 
case. We state and prove the main result. Some of the following lemmas are stated without 
proof; the reader interested in the details is referred to the references. As in the Riemannian 
case the main idea is to construct appropriate bounds for the fundamental amplitudes. 

Lemma 9.3. Ji.9| \2f^ For any subset of k elements p\. . . p K out of the corresponding four 
representations appearing in A e (pi, . . . ,Pa), the following bounds hold: 



\A e (pi, . . . ,p 4 )| < — ^5-, where a 

Up 



C K , f 1 for k < 3 

I for k = 4 



for some positive constant C K . 

Lemma 9.4. \12U^ A e (pi, . . . , p 4 ) and A v (p 1 , . . . , p w ) are bounded by a constant indepen- 
dent of the p 's. 

Lemma 9.5. Wty For any subset of k elements p\. . . p K out of the corresponding ten 
representations appearing in A v (pi, . . . , p 10 ) the following bounds hold: 

K K 

\A v (p!, . . . , pio) I < 



Up, 

i=l 



for some positive constant K K . 



Theorem 9.2. IiM WJJ^ The state sum A(J~), A15ty) . converges for any non degenerate 
2-complex J . 

Proof. We divide each integration region M. + into the intervals [0, 1), and [1, oo) so that 
the multiple integral decomposes into a finite sum of integrations of the following types: 

i. All the integrations are in the range [0,1). We denote this term T(F, 0), where F 
is the number of 2-simplexes in the triangulation. This term in the sum is finite by 
Lemma T 



ii. All the integrations are in the range [l,oo). This term T(0, F) is also finite since, 
using Lemmas 19.31 and 19.51 for n = 4, and k = 10 respectively, we have 



T(0,F)<H J dp f p 2 p n «-^< J dp, pf° 
J Pf =i \pf =1 



< oo. 
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iii. m integrations in [0, 1), and F—m in [1, oo). In this case T(m, F—m) can be bounded 
using Lemmas ()9.3|) and ()9.5|) as before. The idea is to choose the appropriate 
subset of representations in the bounds (and the corresponding values of k) so that 
only the m — F representations integrated over [1, oo) appear in the corresponding 
denominators. Since this is clearly possible, the T(m, F) terms are all finite. 

We have bounded Z by a finite sum of finite terms which concludes the proof. □ 

The 'scaling' behavior of amplitudes in the model is very similar to that of the Rie- 
mannian version of the previous section. If one makes the substitution (2j + 1) = — ip and 
ip y = ii] y in the Riemannian expression for the kernel ()144jl one obtains the Lorentzian 
expression (|152j) . In addition, with the same substitution the face amplitude (2j + l) 2 in 
(j!46j) takes the form p 2 of the face amplitude in (j!53jl . From this viewpoint the finiteness of 
the Lorentzian model is intimately related to the finiteness of the Riemannian counterpart. 

The non compactness of SL(2, C) introduces potential divergences due to the fact that 
volsi(2,c) — oo. This problem can be dealt with in the GFT formulation by using the 
appropriate gauge fixing conditions as explained in Section 0] If we define the model on a 
simplicial decomposition the same gauge fixing can be implemented. 

9.3.3 A new Lorentzian model 

In the Euclidean case there was only one way of selecting the subgroup U of group elements 
leaving invariant a fixed direction in Euclidean space-time. In the Lorentzian case there are 
two possibilities. The case in which that direction is space-like was treated in the previous 
section. When the direction is time-like the relevant subgroup is U — SU(1, 1) x 

This case is more complicated due to the non-compactness of U. Consequently, one 
has to deal with additional infinite volume factors of the form [/-volume in (jll6|) . An- 
other consequence is that the invariant vectors defined in (j!28|) are now distributional and 
therefore no longer normalizable. All this makes more difficult the convergence analysis 
performed in the previous models and the issue of finiteness remains open. 

On the other hand, the model is very attractive as its state sum representation contains 
simple representations in both the continuous and discrete series. As pointed out in [TU] and 
discussed in the following section one would expect both types of simple representations 
to appear in a model of Lorentzian quantum gravity. 

The difficulties introduced by the non-compactness of U make the calculation of the 
relevant kernels (jlHlj) more involved. No explicit formulas are known and they are defined 
by integral expressions. We will not derive these expressions here. A complete derivation 
can be found in |112j . The idea is to use harmonic analysis on the homogeneous space 
SL(2, C) / SU(1, 1) x Z 2 which can be realized as the one-sheeted hyperboloid y^y u r]^ u = — 1 
where y and — y are identified (imaginary Lobachevskian space, from now on denoted H~). 
The kernels correspond to eigenf unctions of the massless wave equation on that space. 
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The corresponding kernel is given by the following expression: 



c+ 

32 Ti-p-^M 6 ^)] 

+ 6(p)S nAk 5{x v i v )5{y V Ql (154) 

where x, y £ and £ G C + is a normalized future pointing null vector in Minkowski 
space-time. The integration is performed on the unit sphere defined by these vectors with 
the standard invariant measure dm. 

As in the previous cases the expressions for the edge and vertex amplitudes (equa- 
tions ()132j) and ()133j) ) reduce to integrals on the hyperboloid H . The expression for the 
amplitude (J133)) of an arbitrary diagram is |112j 

a ( j )=e / d pfH(p 2 f+ n2 f) n MMApe}) n a^kum), 

hi f e v 

where now there is a summation over the discrete representations n and an multiple integral 
over the continuous representations p. For the discrete representations the triangle area 
spectrum takes the form (|97|) if we define k = j+1/2. The finiteness properties of this model 
have not been studied so far. There is a relative minus sign between the continuous and 
discrete eigenvalues of the area squared operators that has been interpreted as providing a 
notion of micro causality in |112j . It would be interesting to study this in connection with 
the models of Section 1531 



10 Discussion 

We conclude this article with some remarks and a discussion of recent results and future 
perspectives in the subject. As we have devoted the second part of this paper almost 
exclusively to the Barrett-Crane model we start by discussing open issues in this context. 
We go to more general open questions toward the end of this section. 

1. Normalization: The Barrett-Crane model is certainly the most well studied spin foam 
model for 4-dimensional quantum gravity. However, as we pointed out at the end of 
Section the definition of the measure is incomplete since there remains to determine the 
value of lower dimensional simplex amplitudes such as the edge (tetrahedron) amplitude 
A e . Although, the normalization (jl451 I14(jj) of the Barrett-Crane model can be fixed in a 
natural way using the GFT formulation — or using similar (subgroup averaging) arguments 
|1181 11211 1122j — one would really like to understand this prescription in relation to the 
continuum theory. 

In Section we have shown that the Barrett-Crane model can in fact be interpreted 
as a 'lattice' formulation of the gravity sector of the 50(4) Plebanski theory — this conclu- 
sion can be extended to the Lorentzian theory without further complications. This close 
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relationship between the model and a continuum action suggested the possibility of de- 
termining the correct spin foam measure by comparison with the formal continuum path 
phase space path integral measure of Plebanski's theory. If the value of the edge amplitude 
(j!45j) could be found in this way this would provide further evidence for the GFT-gravity 
duality and strengthen the fmiteness results. These considerations apply also for other spin 
foam models in particular to those defined by constraining BF theory. A general prescrip- 
tion for the construction of an anomaly free measure for spin foam models is investigated 
in [H3J. 

2. Numerical results: A very attractive property of the normalization (jl46j) is however the 
fmiteness of transition amplitudes on a fixed discretization. Having a well defined model — 
for both Riemannian and Lorentzian gravity — has opened for the first time the possibility 
for numerical explorations. An important step toward the numerical evaluation of the BC 
model is the development of efficient algorithms for the computation of the lOj'-symbols 
(1231 . Numerical calculations by Baez et al. |93] show that the sum over spin foams in 
our normalization of the Barrett-Crane model converges very fast so that amplitudes are 
dominated by spin foams where most of the faces are labeled by zero spin. The leading 
contribution comes from spin foams made up from isolated 'bubbles of geometry'. This 
problem is however not present in the Lorentzian version (|153|) where simple representations 
are in the continuum part of the spectrum. Baez, Christensen, Halford and Tsang propose 
modifications of the normalization which still yield finite amplitudes but avoid this puzzling 
feature. However, these modifications consist of ad hoc definitions of the value of edge and 
face amplitudes. A clear understanding of this would have to come with the understanding 
of the normalization issues discussed above. 

If we take seriously the normalization (|145J) . the previous result might not be so worri- 
some after all. According to Section ITT^ the asymptotic expression of the vertex amplitude 
shows that for large spins it is dominated by degenerate 4-simplexes. These configurations 
are however strongly suppressed by the form of the edge amplitude ()145|) in our normal- 
ization. Assuming that for small spin values the amplitudes behave appropriately, large 
spin behavior might not be physically relevant after all. 

As pointed out by Smolin |124j . one could neglect spin zero representations according 
to the argument that they represent 'zero geometry states' (all geometric operators vanish 
on the spin zero subspace). If one does so, then transition amplitudes are dominated by 
the spin one half representations and the model becomes very similar to a model of dy- 
namical triangulations — in dynamical triangulations one sums over triangulations whose 
1-simplexes have the same definite length, while here the amplitudes are dominated by 
triangulations where triangles have the same fixed area given by (J97|) for j = 1/2. However 
one should also point out an important difference: in our case, configurations are weighted 
with complex amplitudes exp(iS) which no (a priori) connection to the Euclidean ampli- 
tude exp(— S) of dynamical triangulations. In the Lorentzian models the analogy simply 
disappears as representation labels are continuous. 

3. The connection with the canonical picture: We began this article motivating the spin 
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foam approach from the perspective of loop quantum gravity. Can we establish a rigorous 
connection between the spin foam models presented here and the canonical formulation? 
In all the models (for Riemannian gravity) introduced in Section with the exception 
of the Barrett-Crane model, this relationship is manifest at the kinematical level as the 
boundary states are given by S77(2)-spin networks. This kinematic connection exists by 
construction in these cases. 

What happens in the case of the Barrett-Crane model? The Riemannian BC model is 
defined in terms of Spin(4), so naturally one could expect the boundary data to be given 
by Spin{A) spin networks (labeled by two half-integers (je,j r ))- The simplicity constraints 
impose ji = j r = j which could be interpreted as the existence of an underlying SU(2) 
connection. However, it can be shown that the boundary data of the model cannot be 
interpreted as given by an SU{2) connection (recall the discussion of Section 19.1. 4(1 . Sim- 
ilarly, in the Lorentzian Barrett-Crane model boundary states could be naively related to 
SL(2, C) spin networks. However similar considerations as in the Riemannian sector show 
that the simplicity constraints (applied a la Barrett-Crane) reduce the boundary states in a 
way that cannot be interpreted as a boundary connection. Alexandrov |125l H26 has been 
studying the possibility of defining a SL(2, C) connection canonical formulation of general 
relativity. The second class constraints of the theory — related to the simplicity constraints 
of Plebanski's formulation — are dealt with using the Dirac bracket in terms of which the 
original connection becomes non commutative. In this sense there is no genuine or stan- 
dard spin network representation of states. It seems that such approach could lead to the 
understanding of the relationship of the Barrett-Crane model and a canonical formulation 
of gravity. The are in fact some indications of this |127t [128]. 

Another difficulty in making contact with LQG is the fact that operators associated to 
triangle-areas have continuous spectrum in the Lorentzian models which is in conflict with 
the canonical result (|7j). From this perspective one is compelled to study the model we 
have described in Section 19. 3. 31 for which part of the area spectrum coincides with (J7|). 

In the canonical framework the compact SU(2) formulation is achieved by means of the 
introduction of the Imirzi parameter i which defines a one parameter family of unitarily 
inequivalent representations of quantum gravity. In particular, geometric operators are 
modulated by t (see (J7J) for example). Can the spin foam approach say something about 
i? There is no conclusive answer to this question so far. The spin foam quantization 
of generalizations of Plebanski's action including the i ambiguity j!29j are analyzed in 

rnnirm] . 

So far we have only compared spin foams with loop quantum gravity at the kinematical 
level. However, one of the main motivations of the approach was the expected simplifi- 
cation in addressing dynamical questions. As we mentioned before, dynamics is not well 
understood in the canonical formulation. This is mostly due to the difficulties associated to 
regularization ambiguities in the construction of the scalar constraint, the understanding 
of its space of solutions and the the definition of the physical scalar product. As we have 
seen all these issues can be naturally addressed in the spin foam formulation. In turn, it 
would be nice to understand whether the approach can help set a guiding principle that 
would allow for the selection of a preferred regularization of the scalar constraint. Under- 
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standing of this question would strengthen the relationship between the two formulations. 
An explicit analysis of the relation between the covariant (spin foam) and canonical formu- 
lation of BF theory is presented in |132j . The reconstruction of the corresponding scalar 
constraint operator is studied and carried out explicitly in simple models. 

4. Alternatives to the Barrett-Crane model: The fact that degenerate configurations dom- 
inate the asymptotics of the Barrett-Crane vertex amplitude [109J can be interpreted as a 
serious problem. A possible solution to the apparent problem is proposed in |llUj . where 
the evaluation of the vertex amplitude is modified in a simple way to avoid degenerate 
contributions. The potential dominance of degenerate sectors of Plebanski's real formu- 
lation over the gravity sector in a formal path integral has been emphasized on general 
grounds by Reisenberger in |133| . There is so far no rigorous result linking the previous 
two. The spin foam quantization of degenerate sectors (e = 0) of Section 17. lOl is an attempt 
to investigate this problem. A puzzling feature of the degenerate model is that, even when 
the vertex and face amplitudes are different, it fully contains the Barrett-Crane 4-simplex 
configurations. Is this an indication of the apparent problem raised by Reisenberger? This 
motivates the question of whether e in (J72j) is well defined in the Barrett-Crane model, and 
if so, whether it differs from zero. The vanishing of e would appear as a serious obstacle 
to reproducing the gravity sector of Plebanski's theory in the low energy limit. 

5. Coupling with matter. The inclusion of matter into the spin foam formalism is of 
clear importance. Here we review recent results obtained in this direction. A natural 
generalization would be gravity with cosmological constant. The definition of the analog 
of the Barrett-Crane model with cosmological constant has been explored by Roche and 
Noui in [134] . The presence of the cosmological 'matter' is modeled using the quantum 
deformation U q (SL(2,C)) whose representation theory is well understood |135| and the 
deformation parameter is determined by the value of the cosmological constant. In this 
modification the homogeneous space H + is replaced by the quantum deformation of the 
hyperboloid defined in terms of a pile of fuzzy spheres. In a recent paper Oriti and Pfeiffer 
[136 proposed a model which couples the Riemannian Barrett-Crane model with Yang- 
Mills theory. Their construction is likely to be generalizable to other models. In the 
context of the GFT formulation of spin foam models Mikovic proposes a way to include 
matter by adding certain matter fields into the GFT-action. In |137[ I138j he generalizes the 
GFT construction described in Section |S] by allowing for the inclusion of spinnor fields in 
finite dimensional representations of 5*0(2) and 50(3) representing matter corresponding 
to fermions and gauge fields. States in the theory are given by spin networks with open 
links; this is consistent with results obtained in the canonical approach |139t [7] . 

More radical and very appealing possibilities are suggested by Crane. In |142j he pro- 
poses a topological QFT of the type of BF theory as fundamental theory. The gravitational 
degrees of freedom are represented by subset of (constrained) representations while those 
of matter are encoded in the remaining ones. The full spin foam model is topological (no 
local degrees of freedom). In order to recover the low energy world (with local excitations) 
the author appeals to certain "symmetry breaking" of topological invariance. Crane's other 
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proposal consists of interpreting topological (conical) singularities naturally arising in the 
structure of the Feynman diagrams of the GFT theory as representing matter degrees of 
freedom |143j . These possibilities are very attractive for their purely geometric character. 

6. Discretization dependence: In Section|3]we have seen that the discretization dependence 
is trivial in three dimensions. In particular this has been nicely formalized in the definition 
of continuum spin foams by Zapata [70] . Some definition of the refinement limit — discussed 
in Section 16.21 — should be investigated for the models which are not trivially discretization 
independent. The model defined in Section IV. 101 (corresponding to the degenerate sectors 
of Plebanski's theory) might be an interesting candidate. The model is not topological 
and has a clear connection to a continuum action. It is somehow between the theory we 
want to define and the simpler theories we understand well but that do not have local 
excitations (such as BF theory and 3-dimensional gravity). From this viewpoint we believe 
that it might be useful to explore its behavior under refinement as a 'toy model'. The 
simplicity of the model might even allow for analytic computations. Incidentally other 
open issues in the spin foam approach — such as the problem of the continuum limit, gauge 
and anomaly questions and the construction of the generalized projector into the physical 
Hilbert space — should be investigated in this model. 

Another proposal for a discretization independent formulation is given by the GFT 
formulation. The GFT formulation is very attractive since it provides a discretization 
independent formulation of spin foam models from the outset. Also it has been very 
useful for the definition of the Lorentzian models of Sections 19.21 and 19.31 as a device for 
formal manipulations. However, the mathematical consistency of this definition depends 
on whether one can make sense of the expansion in A of equation (jl05|) . Suggestions on 
how the series could be summable by complexification of the coupling A can be found in 
|16j . A beautiful example showing that the A-series can be summed in certain cases has 
been constructed by Freidel and Louapre [HH] for certain 3-dimensional GFT's. 

In the standard background dependent realm, the problem of removing the regulator 
introduced by the discretization is referred to as the continuum limit. Notice that the 
nature of the regulator in spin foams is totally different. The discretization introduces a 
cut-off in the allowed configurations, but the discrete and combinatorial nature of the latter 
is expected to be preserved when the regulator is removed. When the regulator is removed 
we expect to obtain a definition of the generalized projection operator which allows for 
the computation of physical transition amplitudes between spin-network states. On the 
other hand, physical states are expected to lie outside the kinematical Hilbert space and 
so they could be given by elements in the dual that resemble (in some sense) continuous 
configurations. In order to avoid confusion about the matter we have consistently referred 
to the notion of removal of the discretization dependence instead of that of the continuum 
limit. 

The approach of dynamical triangulations is a background independent formulation 
of quantum gravity in which the removal of the regulator is similar in spirit to both the 
problem in the spin foam context and that of (the background dependent) lattice gauge 
theory. Here we present an account of the recent results and point out the conceptual dif- 
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ferences with the spin foam approach. We hope that, despite the differences, these results 
might be helpful in developing useful techniques in the spin foam approach. In dynamical 
triangulations the nature of the regulator is different from that of spin foam models. In 
this approach, (diff-equivalent classes of) smooth metric configurations are approximated 
by space-time triangulations where 1-simplexes have the same fixed proper length £. The 
smaller the length scale the better the approximation; therefore, the proper length rep- 
resents the regulator in the theory to be removed in a certain i — > limit. The phrase 
continuum limit certainly has a clear-cut meaning here. The continuum limit in dynam- 
ical triangulations has been extensively studied. As we previously mentioned, dynamical 
triangulations is a definition of quantum gravity based on the Euclidean path integral, i.e., 
configurations are weighted with real amplitudes exp(— S). In order to recover Lorentzian 
quantum gravity one has to define what is meant by a Wick rotation in the background 
independent framework. Once this is done, the continuum limit can be studied using 
standard techniques of statistical mechanics. Results indicate that in Euclidean dynamical 
triangulations there is no continuum limit |144j . The path integral is dominated by singular 
Euclidean geometries. The situation improves in the Lorentzian models. This is because 
of the restrictions imposed by the notion of Wick rotation defined. These restrictions are 
such that in 1+1 gravity the model is exactly soluble and possesses a continuum limit. 
There is some numerical evidence that could be the case in higher dimensions |145[ I146j . 

7. Contact with the low energy world: In this article we have concentrated mostly on 
mathematical and conceptual issues which are important in the construction of a consistent 
spin foam model for gravity. However, one of the most pressing questions is whether the 
spin foam models can reproduce the low energy world of general relativity. Although this is 
certainly one of the most important questions, there is unfortunately no conclusive evidence 
of this at present and most of the work lies ahead of us. 

Some of the models presented in Section satisfy some notion of 'naive' classical-limit 
in the sense that they are derived from a continuous action of general relativity. Other 
models are defined by postulating the fundamental dynamics and using the kinematical 
structure discovered by the canonical formulation. If any of these models are to be taken 
as strong candidates for theories of quantum gravity they must be able to reproduce the 
physics of gravity at low energy and predict the corresponding semi classical corrections. In 
this respect, it is still an unsettled issue whether the differences of the various models at the 
fundamental level lead to different theories or rather are to be expected to disappear in the 
continuum limit. After all, the latter happens in lattice gauge theory in the renormalization 
process. Even when progress has been made in this direction, there is however no clear-cut 
formulation of the analog of renormalization theory for spin foams. The reason for this is 
the difficulty of applying standard techniques in the background independent context. In 
Section we mentioned some new ideas and attempts to tackle this problem. This is an 
extremely important problem where new ideas will have to play an important role. 

8. The generalized projection vs. the Feynman proper-time propagator. Throughout this 
work we have advocated the viewpoint of the path integral as a device for the computation 
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of the generalized projection operator into the physical Hilbert space of quantum gravity. 
We claimed that this picture is forced upon us by the general covariance of general relativity. 
Although we have kept the term transition amplitude for notational convenience there is 
no notion of time involved in our definitions. 

In other path integral approaches one aims at the construction of the notion of tran- 
sition amplitudes between 3-geometries at a definite time. Such transition amplitudes are 
governed by the so-called proper time evolution operator or proper-time propagator. 

This possibility has also been explored in the context of LQG and is the perspective 
in which the causal models of Section 15.51 are defined. There are also modifications of 
the Barrett-Crane model which attempt to define such proper time evolution operator 
|147| [148 . The main difficulty in such approach is how to give meaning to the notion of 
'time' in the general covariant context. 



The suggestion that quantum gravity should be described in terms of discrete combi- 
natorial structures can be traced all the way back to Einstein |149j . In this review we have 
described an approach inspired by the non-perturbative canonical quantization of general 
relativity in which these suggestions are concretely realized. We have shown that a great 
deal of progress has been achieved in understanding the conceptual issues involved. The 
fundamental motivation was the construction of a theory that would provide a device to 
analyze dynamics in quantum gravity. The conceptual setting is clear: spin foams should 
provide the definition of the physical scalar product and hence the physical Hilbert space 
of quantum gravity which encodes all the information about quantum dynamics. We have 
seen that these models can be obtained as lattice quantizations of general relativity in 
appropriate variables. The 'lattice' action being naively related to the continuum general 
relativity action (or some equivalent classical formulation) in the sense as Wilson's action 
for lattice gauge theory. Some of these models have remarkable finiteness properties both 
in the Riemannian and Lorentzian sectors. We have also shown how the basic structure of 
spin foams arises from the canonical formulation of loop quantum gravity. In addition we 
described models which are not related to a classical action and are constructed from the 
basic properties of spin networks plus simple causality requirements. For all these models 
the common structure arising is given by spin foams: colored 2-complexes where the geo- 
metric degrees of freedom are encoded in a fully combinatorial manner. Spin foam models 
appear as a beautiful realization of Einstein's idea. There are certainly many difficult open 
questions and we have tried to point out those which we judge the most important ones. 
We hope that new ideas and hard work will continue to contribute to their resolution in 
the near future. 
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